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BLACK HOLE EMERGENCE IN SUPERNOVAE 
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ABSTRACT 

If a black hole formed in a core-collapse supernova is accreting material from the base of the envelope, the 
accretion luminosity could be observable in the supernova light curve. Here we continue the study of matter fall 
back onto a black hole in the wake of a supernova and examine realistic supernovae models which allow for an 
early emergence of the accretion luminosity. Such cases may provide a direct observational identification of the 
black hole formed in the aftermath of the explosion. Our approach combines analytic estimates and fully relativis- 
tic, radiation-hydrodynamic numerical computations. We employ a numerical hydrodynamical scaling technique 
to accommodate the diverse range of dynamical time scales in a single simulation. We find that while in typical 
Type II supernovae heating by radioactive decays dominates the late-time light curve, low-energy explosions of 
more massive stars should provide an important exception where the accretion luminosity will emerge while it is 
still relatively large. Our main focus is on the only current candidate for such an observation, the very unusual 
SN1997D. Due to the low energy of the explosion and the very small (2x 10~ 3 M Q ) inferred mass of 56 Co in 
the ejected envelope, we find that accretion should become the dominant source of its luminosity during the year 
2000. The total luminosity at emergence is expected to lie in the range 0.5- 3 x 10 36 ergs s" 1 , potentially detectable 
with HST We also discuss the more favorable case of explosions which eject negligible amounts of radioactive 
isotopes and find that the black hole is likely to emerge a few tens of days after the explosion, with a luminosity 
of- 10 37 ' 



ergs s . 

Subject headings: accretion, accretion disks — black holes 
supernovae: individual (SN1997D) 



methods: numerical — supernovae: general - 



1. INTRODUCTION 

Core-collapse supernovae mark the death of a massive star 
and the birth of a compact remnant. Theory suggests that the 
remnant can be either a neutron star or a black hole, depending 
on the character of the progenitor and the details of the explo- 
sion. Radio pulsar emission has allowed to compile a variety 
of observational evidence associating neutron stars with sites 
of known supernovae, but similar evidence for a black hole - 
supernova connection is still mostly unavailable. Interestingly, 
indirect evidence that the black hole candidate in the X-ray bi- 
nary system GRO J1655-40 was fo rmed in a supernova e xplo- 
sion has recently been reported by Israelian et al. (1999). The 



inference is based on detection of high abundances of nitrogen 
and oxygen on the surface of the companion, which has too low 
a mass to have produced them by thermonuclear burning, hence 
indicating that they were deposited there by the supernova that 
created the black hole. 

The formation of a compact object in a supernova can be in- 
ferred directly if the object causes an observable effect on the 
total luminosity which follows the explosion - i.e., the light 
curve. In particular, some material from the base of the ex- 
panding envelope may remain bound to the compact object and 
continuou sly fall back onto it, thus generating an accretion lu- 
minosity. Zampieri et al. (1998a) have recently performed a 



self-consistent investigation of the accretion luminosity gen- 
erated by spherically symmetric "fallback" of matter onto a 
black hole in the wake of a supernova. With the aid of a 
general-relativistic, radiation hydrodynamic Lagrangian code 
they showed that, at late times, the fallback evolves quasi- 
stationarily with the accretion lumin osity obeying the analytic 
expression found by Blondin (1986) for steady state, spherical 
hypercritical accretion 4 . The accretion luminosity declines sec- 
ularly with a time dependence of f" 25 / 18 . This decay is driven by 
the decrease in the accretion rate due to the continuous expan- 
sion of the supernova envelope (Colpi, Shapiro & Wasserman 
1996), 

Such an accretion luminosity will produce a distinct signa- 
ture on the total light curve if and when it becomes compa- 
rable to the output of other power sources, namely the initial 
internal energy of the envelope and decays of radioactive iso- 
topes synthesized in the explosion. It is well established that 
the efficiency of converting spherical accretion onto a black 
hole into radiative luminosity may be quite low (Shapiro 1973; 
Blondin 1986; Wandel, Yahil & Milgrom 1986; Park 1990; 
Nobili, Turolla & Zampieri 1991; Zampieri, Miller & Turolla 
1996), and, in general, the accretion luminosity will be unde- 
tectable in comparison with the luminosity generated by ra- 
dioactive decays. This is particularly important when compared 
to the case of accretion onto a neutron star, which is expected 
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to be far more efficient in generating an accretion luminosity. 
For the last decade this distin ction has been s pecifically ap- 
plied to the cas e of SN1987A: |Chevalier (1989b and H ouck & 
Chevalier (1991) found that accretion onto a newly formed neu- 
tron star in SN1987A would generate an accretion luminosity 
at the Eddington-rate (;> 10 38 ergs s" 1 ) within a few months af- 
ter the explosion, and the luminosity would persist at this rate 
for several years 5 . The proximity of SN1987A has allowed ob- 
servations to follow the light curve for over 3000 days (to a 
present luminosity of ~ 10 36 ergs s" 1 ), and the light curve has 
been completely consistent with heating by decays of 56 Co , 
57 Co and 44 Ti . The absence of any evidence for accretion lu- 
minosity in the case of SN1987A has been used to argue that the 
compact remnant in SN1987A is a low mass black hole, imply- 
ing si gnificant consequences regarding the nuclear equation of 
state (Bethe & Brown, 1995). It has also been argued, however, 
that this absence could still be consistent with a weakly magne- 
tized neutron star, if the accretion flow has become dynamically 
unstable (Fryer, Colgate & Pinto 1999). 

A positive indication of the presence of a black hole would be 
an actual detection of the accretion luminosity. The power-law 
time dependence would make the accretion luminosity easily 
discernible from luminosity generated by radioactive heating, 
which decreases exponentially with time. Furthermore, the dif- 
ference in the time dependence implies that the accretion lumi- 
nosity must eventually become dominant over radioactive de- 
cays, so in principle, the black hole always "emerges" in su- 
pernovae light curves. However, t he potential for observ ation 
is usually quite slim: for example, Zampieri et al. (1998a) esti- 
mated that if a low mass black hole was formed in SN1987A, 
the accretion luminosity will not become dominant over that of 
heating by 44 Ti decays until ~ 900 years after the explosion, 
when the luminosity will have dropped to a mere 10 32 ergs s" 1 . 
The abundance of radioactive elements and the amount of fall- 
back in SN1987A are rather typical of core collapse supernovae 
of low to intermediate mass progenitors (12-25 M©) so that, if 
a black hole is formed, the "emergence" of the accretion lumi- 
nosity occurs only at such late times and low luminosities that 
it is of no observational consequence. 

An important exception may exist, however, in the explo- 
sions of more massive stars. In their survey of supernov a ex- 
plosions and nucleosynthesis, Woosley & Weaver (1995) find 
that more massive stars tend to produce more massive compact 
remnants: the original iron core is larger, leading to a more mas- 
sive proto-neutron star and to larger amounts of fallback (see 
also Fryer 1999). Correspondingly, larger mass progenitors are 
more likely to produce black holes (as the remnant mass would 
be larger than the maximum mass of stable neutron stars), and 
to have larger amounts of material that remain bound to the 
black hole as a reservoir for late time fallback. Just as impor- 
tant, the final abundance of the key radioactive isotopes ( 56 Ni , 
57 Ni and 44 Ti ) in the expanding envelope is quite sensitive to 
the location of the cutoff that defines the material which settles 
on the black hole during the early stages (while the explosion is 
still proceeding). For the more massive progenitors the cutoff is 
further out, and since the radioactive elements are synthesized 
in the innermost layers above the proto-neutron star, the total 
mass of these elements in the ejected envelopes of higher mass 
stars can be significantly smaller than in "standard" supernovae. 
In fact, progenitor of masses in the range 30-40M Q , are ex- 
pected to be practically free of radioactive elements that can 



power the bolometric light curve ( Woosley & Weaver, 1995 ). 
For such stars, black hole accretion luminosity should emerge 
immediately after the decay of the recombination peak (tens of 
days after the explosions), while it is still relatively powerful, 
and hence potentially observable. 

Naturally, the more massive stars are a small minority in the 
stellar population, as are their explosions among core collapse 
super novae. However, one recen t supernova, SN1997D in NGC 
1536 ( leMello & Benetti, 1997), may present a border-line can- 
didate for a direct detection of black hole emergence. Its ob- 
served light curve indicates that that the total mass of 56 Co (the 
daughter of 56 Ni ) in the envelope is only ~ 2 x 10~ 3 M@, a fac- 
tor of ~ 35 l ower than SN1987A . By analyzing the light curve 
and spectra, [Turatto et al. (1998| ) suggested that SN1997D was 
a low energy (E lot ~ 4x 10 3U ergs) explosion of a 26 M Q star. 
The expected mass of the remnant is about 3 M Q (formed by a 
1 .2 Mq early fallback on a 1.8 M Q collapsed core), hence most 
likely a black hole. Based on the inferred average properties of 
the material available for late time fallback and on the late time 
behavi or of the accretion luminosity found by Z ampieri et al. 
(1998a), Zampieri, Shapiro & Colpi (1998b) estimated that ac- 
cretion luminosity will emerge in SN1997D in the optical band 
within only a few years after the explosion. 

In this work we revisit the issue of black hole emergence 
in supernova due to spherically-symmetric late-time accretion. 
In particular, we investigate a more robust estimate for our 
"proto-type" case, SN1997D. We furnish revised analytic es- 
timates and a detailed numerical investigation of the radiation- 
hydrodynamic evolution of the ejecta, incorporating the effects 
of a realistic chemical composition, opacities and radioactive 
heating on the accretion history and luminosity. The numer- 
ical study is based on an im proved version of the ra diation- 
hydrodynamic code used by Zampieri et al. (1998a), where 
the main modifications are the inclusion of realistic envelope 
composition and opacities (rather than pure hydrogen) and of 
heating by radioactive decays. We find that a 3 M Q black hole 
should emerge in the SN1997D light curve at about 1000-1500 
days after the explosion - i.e., during late 1999 and 2000 - 
when the total bolometric luminosity is ~ 0.5 - 3 x 10 36 ergs s" 1 , 
and is marginally detectable with HST 

In the course of our computational study we were compelled 
to introduce a rescaling scheme in order to numerically follow 
the accretion history and light curve for several years. The 
huge dynamical range of pertinent time scales renders a full 
scale time-integration numerically impractical. Our solution 
has been to study a scaled model which evolves more quickly 
and allows a simple relation between the properties and his- 
tory of the rescaled model and the true one. We note that this 
scheme, which is presented in detail in the appendix, is poten- 
tially suitable for other problems which involve time-dependent 
radiation hydrodynamics spanning a large dynamic range. 

We begin our discussion in § || with a summary of the char- 
acteristic features of supernova light curve evolution when ac- 
cretion onto a central black hole is occurring. Quantitative re- 
lations between the light curve magnitude and the relevant time 
scales are al so reviewed. The princ ipal features of the numer- 
ical code of Zampieri et al. (1998a) and the modifications in- 
corporated for this work are described in § 0. The light curve 
of SN1997D is examined in detail in § |], where we present 
analytic and numerical results, including our estimates regard- 
ing black hole emergence. Prospects of observing black hole 



5 It is noteworthy that magnetic dipole emission of a Crab-like pulsar would deposit energy in the envelope at a similar rate ( [Woosley et al., 1989 ). 
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emergence in other supernovae are discussed in § | 
discussion and conclusions are offered in § 0. 



, and further 



2. THE SUPERNOVA LIGHT CURVE INCLUDING ACCRETION 

The bolometric luminosity that follows the explosion of a 
massive star is powered by the internal energy of the expanding 
envelope, emitted as thermal photons. Under "bolometric" we 
do not include the hard X-ray and 7-ray photons that are emit- 
ted in radioactive decays and escape from the envelope before 
thermalization. In most core collapse supernovae, luminosity in 
the early (up to tens of days) light curve is powered by energy 
deposited in the envelope by the supernova shock. After this en- 
ergy is depleted, only ongoing sources of energy can continue 
to power the later part of the light curve, often referred to as a 
"tail". By definition, radiative emission from late fallback qual- 
ifies as an ongoing source, as do decays of radioactive elements 
in the envelope 6 . 

It is useful to examine the evolution of the light curve 
in terms of the key characteristic time scales and luminosity 
scales, which can be estimated on the basis of the explosion en- 
ergy and the average properties of the progenitor (Arnett 1980, 
1996). When accretion onto the central object caused by late- 
time fallback is taken into account, additional time scales and a 



luminosity scale must also be included ( Zampieri et al., 1998a). 



We briefly review these basic estimates below. They will serve 
also as the basis for establishing the rescaling scheme presented 
in the appendix. 

2.1. Expansion 

For an envelope with initial radius R$ and outer velocity Vb 
we can define an initial expansion time scale, to, as 



Vo 



(1) 



Assuming that after the passage of the shock the envelope set- 
tles into free streaming and homologous expansion (as is typi- 
cal in strong spherical shocks), the initial velocity of a shell at 
radius r in the expanding envelope is simply vo(r) = r/to. For 
free streaming the radius of the shell will follow r~?, and since 
the expansion is nearly adiabatic, the temperature T and density 
p of the radiation-pressure dominated gas decline according to 
T ~? and p~f~ 3 . These quantities determine the leakage rate 
of thermal photons from the envelope material, and hence gov- 
ern the history of the early light curve. 



a homogeneous expandi ng envelope t his diffusion luminosity 
can be approximated by (A.rnett, 1996) 



Ldiff(t) = Ldiff,oe 



-(tot+t 2 /2)/(t t d , fffi ) 



where 



Ldiff,o : 



(f7r/3ac)7oVo 



tdifffi : 



(3c 



(2) 



(3) 



with j3 a numerical factor which depends on the temperature 
profile. The value (3 = 13.8 (exact for the "radiative-zero" tem- 
perature profile, ro(r) = 7()(sin(7r^-)/(7r^)) 1 / 4 ) is usually a rea- 
sonable estimate (Arnett, 1980). The time dependence of equa- 
tion (Q) includes the expansion time fo and the initial diffusion 
time tdif/fl, which is the characteristic time scale for radiation 
to cross the envelope through diffusion at the onset of expan- 
sion. The characteristic time for radiation to diffuse out of the 



expanding envelope is given by (Zampieri et al., 1998a) 
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The diffusion approximation breaks down when expansion 
makes the envelope transparent to its thermal photons, after a 
typical time 
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which is significantly longer than tuff. However, the onset of 
recombination in the envelope usually occurs early enough so 
that the actual value of t trans is irrelevant. 

2.3. Recombination 

As the temperature of the envelope is gradually degraded by 
expansion, it eventually reaches the characteristic recombina- 
tion temperature of the envelope material. In a hydrogen-rich 
envelope this temperature can be estimated by T reCi H (approx- 
imately 10 4 °K), so that for adiabatic expansion we find a re- 
combination time scale 



tree - to- 



(6) 



2.2. Diffusion 

During the time that the initial internal energy is the domi- 
nant source for radiative luminosity, the evolution of the light 
curve is governed by the properties of the optically-thick, 
hydrogen-rich outer layer, that holds most of this energy. We 
denote the mass, initial average density and initial average tem- 
perature of this layer as Mg, po : H and Tqh, respectively. 

The initial post-shock temperature of the envelope is gener- 
ally high enough (7b.// > 10 6 °K) so that the envelope mate- 
rial is highly ionized, and the dominant source of opacity, k, 
is Thomson scattering. Correspondingly, the thermal photons 
are mostly trapped in the envelope, and a radiative luminosity 
occurs through diffusion of thermal energy to the surface. For 



The main impact of recombination in the context of the light 
curve is that it imposes a sharp decrease in the opacity to the en- 
velope's own thermal photons: as the envelope material recom- 
bines, it rapidly transforms from opaque to transparent. This 
process is rapid enough so that the photons cannot readjust to 
the global thermal profile of the envelope; rather, a recombina- 
tion "front" sweeps inward through the envelope (but outward 
in the laboratory frame), essentially liberating all the thermal 
energy which is still stored in the envelope at the time t rec . As 
the envelope material recombines, it rapidly transforms from 
opaque to transparent. In many Type II supernovae the time 
until adiabatic expansion cools the envelope to the recombi- 
nation temperature is usually smaller or equal to the diffusion 

s Ybung & Branch (19891 pointed out that if the energy generation from radioactive decays is sufficiently high and the progenitor radius is sufficiently small, 
radioactive decays can become the dominant source of power for the light curve even very early after the explosion. They suggest that this is the case in Type II-linear 
supernovae. 
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time scale and significantly shorter than the transparency time 
scale, so a significant fraction of the initial thermal energy is 
still available at the time of recombination. This energy is then 
released over the time required for the hydrogen envelope to 
recombine, generating what is usually an observable peak in 
the light curve (SN1987A providing a proto-typical example). 
Roughly, the average luminosity during the recombination peak 
can be estimated as 



Ethermal(l rec) 47T 



t^2 7^2 
n 3 1 O.H 1 rec.H 
= KnCL : — 

3 to 



(7) 



assuming that energy losses due to photon diffusion have been 
negligible, which is a good approximation if f ree ,fl < Uiff- In 
equation (0) a is the radiation constant. Note that the nonlinear 
nature of the physics governing the propagation of the recom- 
bination front through the envelope limits the applicability of 
using the average quantities of the envelope to fully describe 
the light curve during the recombination phase (see more de- 
tailed estimates in Nomoto et al. 1994, Arnett 1996). 



2.4. Radioactive Heating 

A complete analysis of a realistic light curve must also in- 
clude the effects of heating due to decays of radioactive iso- 
topes synthesized in the explosion. After recombination, we 
can estimate that the envelope is practically transparent to its 
thermal photons, so that the luminosity due to radioactive heat- 
ing is roughly equal to the instantaneous energy deposition rate 
due to the decays, through thermalization of the decay products. 
The rate of energy deposition in the envelope by the decays of 
an initial total mass Mx of a given isotope X can be expressed 
as: 

Q x (t) = M x (e Xn fx(.t) + exs) z~" Tx , (8) 

where tx is the life-time of the element X. The form of 
equation (||) distinguishes between the photons (7-rays) and 
positrons emitted in the decays (with energy rates per unit 
mass of ex. 7 and ex, e + respectively) since 7-rays are not to- 
tally trapped in the envelope, and their contribution to heating 
is modified by a trapping factor, /(f). On average, the hard pho- 
tons lose about half their energy in every scattering, until their 
energy is degraded to tens of ke V. At these energies free-bound 
absorption on heavy elements becomes the dominant source of 
opacity and the photons are thermalized rapidly. Because the 
7-ray opacity rises very rapidly with decreasing energy, we can 
assume to first order that once a hard photon scatters it is ab- 
sorbed. The trapping factor can be simply approximated as 
(|Woosley et al„ 1989|) 



fx(t)= 1-exp 



-K 



(9) 



where <l>o(f / to) 2 and k 7 .x are the 7-ray column depth and opac- 
ity for the typical photons emitted in the decays of element X. 
Since the inner (more dense) layers are the dominant contrib- 
utors to the optical depth for 7-rays, the helium-rich layer can 
be used to estimate the total 7-ray optical depth. The time- 
dependent column depth is 
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pdr « p He (t)RHe(t) = PHefl ~ 
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In equations (|9]-|To|) it is implicitly assumed that the expansion 
time fo is universal to the entire envelope; otherwise, it must be 
replaced with the characteristic expansion time of the helium- 
rich layer. Note that equation (|l0|) allows us to express the 
trapping factor in terms of the initial (and presumably known) 
global quantities of the helium-rich layer. The quantity po.Heto 
(which appears here to the power |) is a useful one also in the 
context of accretion (see Chevalier 1989, and in detail below). 

In addition to the characteristic time Tx, another character- 
istic time involving radioactive heating is the 7-transparency 
time, similarly to equation (^|). Since 7-rays carry the bulk of 
the energy emitted in the radioactive decays, it is important to 
compare this time, 
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with the other characteristic times governing the light curve his- 
tory. At times t ;> t, rm - 7 ,x some of the hard photons emitted 
in the decays of element X escape before thermalizing and do 
not contribute to the bolometric light curve. We note again, 
that equation ( |TTj ) assumes, as before, that the bulk of the enve- 
lope's opacity to hard photons is contributed by the helium-rich 
layer. The relevant radioactive isotopes for powering the light 
curve are 56 Ni and its daughter nucleus 56 Co , 57 Co and ^Ti . 
We adopt the charac teristic paramet ers f or these isotopes from 
Woosley et al. (1989j ), listed in table ^4] Note that in a realistic 
application, especially when the early history of the envelope 
is of interest, the 56 Co mass m ust be adjusted according to its 
production rate in 56 Ni decays ( [Shigeyama & Nomoto, 199C ), 



M 5f , Co (f) = M 5 6 Ni (t = 0) [exp(-f/T5 f , Co )-exp(-f /T56 Ni )] . (12) 

TABLE 1 

The relevant radioactive isotopes for powering a supernova light curve: 
Energy emission rates in 7— rays and positrons, life-times and 
effective 7— ray opacities at the characteristic photon energies 8 



isotope 


e 7 


£ e + 


T 


^/Yj b 




(ergss ') 


(ergs s~') 


(days) 


(cm 2 /gm) 


56 Ni 


3.90 x 10 10 





8.8 


0.06 


56 Co 


6.40 x 10 9 


2.24 x 10 s 


111.3 


0.06 


57 Co 


6.81 x 10 8 





391.0 


0.144 


44 Ti 


2.06 x 10 s 


6.536 x 10 7 


3.28 x 10 4 


0.073 



"adapted from Woosley et al. 1989. 
^Yj = the total fraction (bound and free) of electrons per nucleon 

2.5. Accretion 

Late-time fallback onto the black hole can be characterized 
in terms of the ratio of the initial expansion time scale to and 
the initial accretion time scale, 



tacc.O — 



GM BH 
s 3 o(He) 



(13) 



where Mbh is the black hole mass and c s .o(He) is the initial 
sound speed in the helium rich layer, which we assume is the 
source of material for late-time accretion. The hierarchy of 
these two time scales determines the hydrodynamic evolution 



5 



of the accretion flow ( 2olpi et al., 1996 ). If k = f flC r.oAo 3> 1, 
the gas has no time to respond to pressure forces, and accre- 
tion proceeds in a dust-like (pressure-free) manner from its on- 
set. Alternatively, if k <C 1 initially, accretion is at first almost 
unaffected by expansion and follows a sequence of Bondi-like 
(Bondi 1952) quasistationary states with a slowly decreasing 
density at large distance. Even in this latter case the flow does 
eventually become dust-like, since expansion continuously de- 
creases the density and pressure in the reservoir of bound mate- 
rial. The transition is expected when the actual expansion and 
accre tion time scales be come comparable, yielding a transition 
time ( |Colpi et al., 1996| ) 



taccX) 



i/(9r-ii) 



£-9(r-i)/(9r-ii) 



(14) 



where V is the adiabatic index of the gas. Thus t, r sw |& 3 f aeCi o 



for T= 4/3 andr„ 



<9\-l/47- 3 /2 



k- i ' 2 t accfi forT = 5/3. 



The hierarchy of time scales reflects the corresponding hier- 
archy of radii, when comparing the marginally bound radius, 
R m b with the accretion radius, R acc . Initially 



Rmbfl = (2GM B Hto) 1 ) Race, = 



GM B 



,(He) 



(15) 



so that if R acc .Q <C R m bfi, the initial accretion is Bondi-like. As 
the e nvelope expands R„ rr ~ t while R„,h ~ t 1 ^ (C olpi et al., 
1996), so that, as in the case of time scales, eventually the ac- 
cretion will become dust-like, even if it had begun as a fluid 
one. The transition time t tr roughly corresponds to the time 
when R mb = R acc . 

If t acc X)/k) 3> 1, pressure forces are never dominant and the 
late time accretion rate from an homologously expanding enve- 
lope can be estimated by the dust-like solution of C olpi et al. 
(1996). They find that the ratio of the dust-accretion rate M to 
the Eddington accretion rate M£dd (see below) is 



M 
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(16) 



by fallback from the expanding envelope closely follows the 
steady-state, spherical hypercritical accretion formula derived 
by Blondin(1986), 
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(17) 



M, 



Edd 



for a given accretion rate M. In equation ( [171 ) n is the opac- 
ity and \i is the mean molecular weight per electron of the ac- 
creting material. The luminosity and accretion rate in equa- 
tion do) are measured in units of the appropriate Edding- 
ton quantities: the Eddington luminosity, L^dd = 4ttGMbh/k = 
1.3xlO 38 (re/O.4cm 2 g _1 ) _1 (M BH /M ) ergss -1 , and the Edding- 
ton accretion rate, M^dd = l^Eddlc 1 . The evolution of the late 
time accret ion luminosity for a du st-like flow can then be ex- 
pressed as (Zampieri et al., 1998b) 
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where 
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OA 



-i/2 / M, 



Mr 



2/3 



5/6 20/9 (19) 



and Aw 1.27 x 10 ergs s (for po.He in gm cm -3 and to in 
seconds) arises from the radiative efficiency of equation (|l7]). It 
is this power-law decay rate in the bolometric luminosity which 
characterizes the presence of a black hole in the aftermath of a 
supernova explosion. 



It is important to notice that equations (18) and (Il9j) cannot 



The resulting luminosity produced by late-time fallback can 
be computed by including radiation in the calculation. A self- 
consistent analysis of the raof/af/on-hydrodynamic evolution of 
the ac cretion flow from an expanding envelope (Z ampieri et al., 
1998a) shows that, at sufficiently late times, the evolution of the 
flow always proceeds as a sequence of quasistationary states. 
This is because the dynamical time scale at the accretion or 
marginally bound radius is much longer than all the relevant 
time scales of the flow and the radiation field in the inner ac- 
creting region. Both pressure and radiation forces are negligi- 
ble at late times, and the flow is indeed dust-like and accretion 
rate declines as f" 5 / 3 . If initially W,oAo ;> 1 and radiation pres- 
sure is always negligible throughout the evolution, M is given 
by equation (|6|) with the initial characteristic parameters of the 
bound material. On the other hand, if t acc o /to < 1 and the ac- 
cretion begins as a fluid flow, the late time accretion will still 
settle on a dust-like solution with M cx t " 5//3 , but with a modified 
(reduced) absolute magnitude. 

As a consequence of the quasi-stationary evolution of the 
flow and the radiation field, the late-time luminosity produced 

7 In the case of a pure dust, the build-up time can be estimated by comparing the 
to the late-time one. We find that the two become comparable at a time of ~ 1.2fo- 



apply at arbitrarily early times even if the initial conditions sat- 
isfy taccf, > h). First, some build up time is naturally required 
for the accretion rate to reach a maximum and start decaying. 
This time should be comparable to the shorter of the two time 
scales to and f acC| o 1 . Furthermore, if during the early evolution- 
ary stages the accretion rate can be so large that the accretion 
luminosity approaches the Eddington limit, radiation pressure 
will modify the flow even if pressure forces were initially neg- 
ligible. On the basis of equation (18), we can define a critical 
time t C rit when the luminosity due to dust-like accretion equals 
the Eddington limit 



fl \ "24/25 , K v 9/25 / MgH 

3 /5 ,8/5 



6/25 



x (20) 



PQ,He 



10 4 gm cm 



(tg.HeY 

V 1 hr/ 

According to equation ( |17| ), the critical rate at 
LEdd is M crlt J> IMq y^ for black holes of 



hrs. 



which L acc ~ 
of several solar 
masses. If the luminosity reaches the Eddington limit, radi- 
ation pressure cannot be ignored in comparison to the grav- 
itational pull of the central object, and so material near the 
marginally bound radius (where a significant fraction of the 
initially bound mass resides) may receive a sufficient impulse 
to become unbound. Hence, we expect that if the build-up 

early time accretion rate, which builds up as M oc t, (Colpi et al. 1996 (eq. [26])) 
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time of the flow satisfies min(fo,f flC c.o) "C f C n>, radiation forces 
will modify the accretion history and limit the accretion rate 
so that M(t J> tent) ~ M cr it by effectively readjusting the val- 
ues of po,He an d f o m the bound region. The total amount of 
bound material must be reduced as well. If, on the other hand, 
min(fo,?flcc,o) 3> tcrir we can expect that radiation pressure will 
be of lesser significance throughout the accretion history and 
that, at late times, the accretion rate will settle on the dust solu- 
tion of equation ( |l6[ ) with the original values of pope and to (if 

ta.o/to ^ !)■ 

2.6. Black Hole Emergence 

We define the black hole "emergence time", tsH , in a super- 
nova light curve to be the time when the accretion luminosity 
is comparable to all the other sources of luminosity combined, 

i-e-, L acc (tBH) = jL ror (tBH)- 

Even a small abundance of radioactive isotopes is sufficient 
to impose a bolometric luminosity that will initially dominate 
accretion luminosity, since the latter is inherently limited by 
the Eddington rate. However, continuous spherical accretion 
must eventually result in emergence, since the accretion lumi- 
nosity (eq. [1180) decreases as a power law, while radioactive 
heating decreases at least exponentially, and even more rapidly 
if 7— ray trapping is incomplete (eq. ||]). Quantitatively, we 
expect that for any specific radioactive isotope, at a time of sev- 
eral times its life-time, tx, the heating rate is declining much 
more rapidly than the accretion luminosity. If the order of mag- 
nitude of the accretion luminosity and initial abundance of the 
dominant radioactive isotope, Mx are known, the time of emer- 
gence of the accretion luminosity is determined mostly by the 
value of Tx, with some finer tuning due to the exact values of 
L aC c(t), Mx, and the time dependence of 7-ray trapping. Note 
that the efficiency of 7-trapping decreases with time as e (_1//f \ 
so the dependence of the time of emergence on Mx and L acc is 
even weaker. Correspondingly, the luminosity at emergence is 

—25/18 

mostly determined by L acc ( several tx) oc L acc ,oT x 

In accordance with § 2.5, it is reasonable to assume that after 
a time not longer than a few x min(fo,f aC e,o) the accretion lu- 
minosity is modulated to L acc < L^dd- Hence, for min(fo,f acCi o) 
of the order of a few days, it is evident that for "typical" core 
collapse supernova, like SN1987A, with an inferred M_w Co ~ 
0.1 M G and assumed Mu Ti ~ 10" 4 M© (Woosley & Weaver, 
1995), the relevant isotope for emergence is 44 Ti , and the 
time scale for emergence is tgH ~ several tu Ti ^hundreds of 
years, at which time the accretion luminosity will drop to a few 
10 32 - 10 33 ergs s" 1 . Note that since about | of the decay energy 
of 44 Ti is emitted in the form of positrons, 7-ray transparency 
will not affect these estimates significantly. On the other hand, 
in a supernova where the abundance of radioactive elements is 
reduced by a factor of ~ 35, as is expected in SN1997D (see be- 
low), 56 Co and 57 Co are also relevant for examining emergence. 
The time of emergence is reduced to hundreds of days, and the 
luminosity at emergence will be a few 10 35 - 10 36 ergs s" 1 . Of 
course, if no radioactive elements are present, as is expected if 
the mass cutoff in the supernova is very far out, emergence can 
occur as soon as the recombination peak clears the envelope, at 
50- 100 days, with L acc > 10 37 ergs s" 1 . 

3. EQUATIONS AND NUMERICAL METHOD: SUMMARY AND 
MODIFICATIONS 

In order to explore the possibility of detecting black hole 
emergence in realistic supernovae in detail, we perform a nu- 



merical study. Time-dependent simulation of the supernova en- 
velope including accretion requires solution of the relativistic 
radiation hydrodynamic equations coupled to the moment equa- 
tions for radiation transport and makes use of a detailed model 
for the ejecta that includes radioactive decay, realistic compo- 
sition and opacities. The equations of relativistic, radiation hy- 
drodynamics for a self-gravitating matter fluid interacting with 
radiation and the general relativistic moment equations for the 
radiation field have been presented in Zampieri et al. (1998a) (to 
which we refer for details). These equations were solved using 
a semi-implicit Lagrangian finite difference scheme in which 
the time step is controlled by the Courant condition and the re- 
quirement that the fractional variation of the variables in one 
time-step be smaller than 10%. 

A fundamental obstacle in tracking the evolution of envelope 
expansion and accretion for several years arises from the ex- 
tremely wide diversity of relevant time scales in the problem. 
These range from the dynamical time in the accreting region 
(milliseconds close to the horizon) to days for the expansion 
time scale, to hundreds of days for the evolution of the light 
curve. 

Since the innermost layers of the envelope have the shortest 
physical time scales, the computational time can be consider- 
ably reduced with a careful choice of the location of the inner 
boundary. During computation, we keep the inner radial bound- 
ary of the integration domain R,„ at hundreds to thousands of 
Schwarzschild radii, far outside the black hole horizon, but still 
but always well within the accreting region. This measure is 
acceptable because the gas in the inner accreting region is in 
local-thermal-equilibrium (LTE) and in near free-fall, thus al- 
lowing a rel iable set of boundary conditions at /?,-„ (Z ampieri 
et al., 1998a). However, during evolution the effective optical 
depth in the accreting region decreases because of the secu- 
lar decrease in density caused by expansion. Thus, in order 
to ensure that the gas is always in LTE at /?,■„, the inner radial 
boundary must be continuously moved inward during the sim- 
ulation. We note that with this treatment, the position of the in- 
ner boundary usually limits the effects of general relativity to be 
quite small. However, we did not alter the general-relativistic 
nature of the code, especially since we cannot determine a priori 
how far out the inner boundary can be set (in particular at late 
times, when it is continuously being moved inward). Further- 
more, while the accretion luminosity is insensitive to placing 
the inner boundary away from the black hole horizon for the 
canonical parameters of SN1997D, in more general scenarios it 
may be essential to apply the exact general relativistic formal- 
ism with an inner boundary in the strong-field domain. 

Even taking advantage of the increase in the physical time 
scales (and hence in the time step) by suitably positioning 
Ri„, the computational time needed to evolve the solution up 
to emergence is excessive. To reach the emergence stage in 
a reasonable computational time, the code employs a "Mul- 
tiple Time-step Proce dure" (MTP), originally developed by 
Zampieri et al. (1998a). With this technique, the integration do- 
main is divided into several sub-grids, and each is evolved on its 
own characteristic time step. Starting from the same initial con- 
ditions and integration domain, the simulation performed using 
this MTP proved to be faster by a factor of 5 - 8 than one where 
no time-step acceleration algorithm was employed. We note 
that the most sensitive aspect of the MTP is maintaining accu- 
rate communication between the subgrids; correspondingly, it 
can be used safely only when there are no sharp features cross- 
ing the envelope (i.e., after recombination). 
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The main focus of the preliminary investigation of Z ampieri 
et al. (1998a) was to study the fundamental features of fall- 
back in a supernova and, in particular, to determine the gross 
properties of the accretion history and luminosity. To this end, 
a simplified pure-hydrogen model of the envelope, generally 
sufficient to produce a typical early light curve of a Type-II 
supernova, was adopted. Here our goal is to examine black 
hole emergence in realistic supernovae, a nd so we expand the 
original survey of Zampieri et al. (1998a ) to take into account 
for variable chemical composition of the expelled envelopes, 
namely radial dependent abundances of various elements. In 
particular, we must account for the presence of radioactive iso- 
topes and the heating of the envelope by their decays. 

3.1. A Variable Chemical Composition 

Inclusion of a variable composition is necessary to reproduce 
realistic superno va light curves, as has been established in stud - 
ies of SN1987A (IWoosley, 19881; Shigeyama & Nomoto, 19901). 



For the purpose of following the radiation hydrodynamic evo- 
lution of the accretion flow onto the black hole, chemical com- 
position and opacity in the accreting region can have important 
quantitative and qualitative effects. First, there is an inherent 
quantitative dependence in the estimate of equation (|l7|), which 
suggests that a pure hydrogen envelope produces an accretion 
luminosity larger by a factor of a few than a helium-metal-rich 
one. Furthermore, it is possible that the presence of metals in 
sufficient quantities may impose a qualitative change in the ac- 
cretion history due to increased opacity where the material is 
partially ionized. Additional effects may arise through the the 
equation of state, but in the case of supersonic accretion we 
expect these to be negligible. 

In this work we study a supernova envelope with a mixtures 
of hydrogen, helium and oxygen, where the relative partial frac- 
tions vary along the radial profile of the envelope. Oxygen is 
expected to be the most abundant heavier element in the ejected 
envelope, and we adopt an approximation that it represents the 
entire metal component in the envelope's final composition (we 
note that over most of the density and temperature ranges of 
interest, an addition of iron at solar abundance would change 
the frequency-integrated opacities by no more than 50%). For 
opacity and ionization fracti ons, we use the TO PS opacities 
available at the LANL server (Magee et al., 1995), with extrap- 



olation to low temperatures based on Alexander & Ferguson 
(1994). The data are stored in the form of tables, with spacing 
of five entries per decade in both temperature and density for a 
given composition. In the course of computation the Rosseland 
and Plank opacities and ionization fractions are calculated with 
linear interpolation in temperature, density and chemical com- 
position. The gas equation of state is approximated as an ideal 
gas of ions and electrons. 

3.2. Radioactive Heating 

For the purpose of computing the contribution of radioac- 
tive decays to the light curve, a consistent treatment of radioac- 
tive heating is required to follow the trapping efficiency, / 7 , 
throughout the envelope. If the rate of radioactive heating is 
large enough it may also affect the hydrodynamic evolution of 
the envelope material (there exists observational evidence for 
such an effect by 56 Ni decay in the early history of SN1987A 
- the "Ni-bubble", Arnett et al. 1989). The role of radioac- 
tive heating may be especially important in the case of inter- 
est here, since most of the radioactive elements are likely to 
lie in the inner part of the envelope, which is the reservoir for 



accretion. The exact distribution will depend on the extent of 
mixing cause d in the process of the explosion ( Shigeyama & 
Nomoto, 1990; pethe, 1990[ ). Significant heating at early times 
may accelerate some of the material which was initially bound 
to the black hole and unbind it, hence imposing a smaller late- 
time accretion rate. On the other hand, at late times, even a 
small amount of heating can maintain a finite ionization frac- 
tion, which in turn will impose a larger photospheric radius, 
hence affecting the luminosity. 

In our numerical study we introduce an effective treatment 
of radioactive elements by including a local heating rate. We 
assume that the abundance of radioactive elements is propor- 
tional to the oxygen fraction ( Pinto & Woosley, 1989 ), with the 
relevant radioactive isotope s lis ted in table |2.4[ We use the ap- 



proximation discussed in § 2.4 that a 7-ray photon is absorbed 
after its first scattering. The trapping factor of the 7-rays emit- 
ted at radius r in the envelope can then be estimated as 



fx, 7 (r)= l-exp<-( 



Kx^pi^dr) 



(21) 



where the quantity in parenthesis is the appropriate 7-optical 
depth at radius r. In our simulation we impose a further ap- 
proximation that all the trapped energy is absorbed locally at 
the point of emission. We note that this simplified treatment of 
radioactive decay energy deposition is satisfactory during the 
first evolutionary stages when the ejecta are still very optically 
thick to 7-rays. Later on, however, the constant opacity/single 
scattering assumption gives only an approximate description of 
the transfer of 7-rays through the envelope. In fact, as the ejecta 
become marginally thick, photons can be scattered or absorbed 
after traveling a significant distance, making non-local effects 
important. An accurate calculation of the energy deposition and 
the 7-ray luminosity would require the solution of the energy- 
dependent transport problem for the 7 rays, which we do not 
attempt. Our effective approach reduces the treatment of ra- 
dioactive heating to an additional local energy source. By de- 
termining the total deposited energy rate 



(22) 



through equations (^J) and (|1]), we modify the energy equation 
(eq.[25] of Zampieri et al. 1998a) by adding the energy input 
from radioactive decay 



e tt + akp(B - wq) + p 



Q, 



(23) 



where P, p and e, are the pressure, mass density and internal 
energy per unit mass of the gas, a is related to the 00 coeffi- 
cient of the spherically symmetric, comoving-frame metric, k p 
is the Plank mean opacity, B = aT A and 4-7rwo is the radiation 
energy density (see Zampieri et al. 1998a for details). All these 
quantities are measured in the comoving frame. 

With the aid of this modified numerical code and our ana- 
lytic estimates we now proceed to study the evolution and light 
curve of realistic supernovae including fallback onto a nascent 
black hole. Our emphasis in this work is the peculiar Type II 
SN1997D, that may provide the first direct observational evi- 
dence of black hole formation in supernovae. This source is 
the only current candidate where an actual observation of the 
light curve at emergence of the black hole may be (marginally) 
possible. 
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4. BLACK HOLE EMERGENCE IN SN 1 997D 

SN19 97D was serendipitously d iscovered on January 14, 
1997 by deMello & Benetti (1997 ) during an observation of 
the parent galaxy NGC 1536. I n their analysis of th e bolo- 
metric light curve and spectra, Turatto et al. (1998 ) found 
that SN1997D was an extremely peculiar Type II supernova. 
Modeling of the early light curve indicates that the super- 
nova was discovered close to maximum light (hence, presum- 
ably, at the recombination peak). Adopting a distance modu- 
lus of (i= 30.64 (a distance of ~ 14 Mpc to NGC 1536), im- 
plied that the maximum absolute magnitude of SN1997D was 
My ^-14.65 (Lmax <; 10 41 ergs s" 1 ), fainter by about 2 mag- 
nitudes than typical Type II supernovae at maximum (P atat et 
al., 1994): SN1997D is possibly the most sub-luminous Type II 
supernova ever observed. 

4. 1 . Fundamental Estimates about the character ofSNl 997D 

The most unique feature of SN1997D in the context of black 
hole emergence is the very low mass of 56 Co deduced from the 
light curve tail (beginning at ~ 30 days after the first observa- 
tion). The magnitude and shape of the tail are consistent with 
a very low abundance of 56 Co , M56f n »2 x 10~ 3 M CA (T uratto 
et al., 1998). This impo rtant result has been recently confirmed 
by Benetti et al. (200(jj ) on the basis of a new set of photometric 
data which extends up to 1.5 years after the explosion. 

turatto et al. (L99is) found that the light curve and spectra 
were best reproduced with 4x 10 50 ergs explosion of a 26 M Q 
main sequence star (see Young et al. 1998 for details). The 
progenitor structure prior to collapse was an 8 M Q helium-rich 
core surrounded by an 18 M hydrogen-rich envelope. The 
remnant formed in the supernova is expected to have a total 
mass of about 3 Mq, composed of a 1.8 M Q collapsed core 
and about 1 .2 M of material which was accreted in early fall- 
back (T. R. Young 1999, private communication). The large 
amount of fallback is a consequence of the low explosion en- 
ergy, and provides a natural explanation for the very low ob- 
served 56 Co mass. Since 3 M Q is almost certainly larger than 
the maximum mass of neutron stars (even a rapidly rotating 
one; Cook, Shapiro & Teukolsky 1994), we can expect that a 
black hole was formed by SN1997D. 

We note that a differen t model of a low mass (~ 8 M Q ) pro- 
genitor was suggested by Chugai & Utrobin (1999 ) to explain 
the early light curve, as well as a low '"Co abundance. In 
such an explosion little fallback would have occurred and the 
remnant would be most likely a neutron star. The late time 
l ight curve shows no evidence of the presence of a neutron star 
(Benetti et al., 2000), which would have been observable if it 



were a Crab-like pulsar, or if it was spherically accreting from 
late-time fallback. Hence a black hole (and therefore a mas- 
sive progenitor) is more consistent with the existing photomet- 
ric data. 

Assuming that the black hole model for SN1997D is cor- 
rect, its remnant includes a black hole surrounded by a low- 
velocity ejecta with a small abundance of radioactive elements. 
If late-time fallback is proceeding in spherical accretion, the 
black hole could emerge rather early, while it is still potentially 
detectabl e. With a preliminary analytic estimate, Z ampieri et 
al. (1998b) found that emergence could occur as early as three 
years after the explosion and that the accretion luminosity at the 
time is roughly in the range 10 35 - 5 x 10 36 ergs s -1 . 

We have reconsidered black hole emergence in SN1997D to 
determine if the accretion tail persists in a realistic model of 



the ejecta that includes a variable composition, realistic opaci- 
ties and radioactive heating and takes into account the radiation 
hydrodynamic evolution of the helium mantle starting from the 
conditions at few tens of hours after the explosion. In the fol- 
lowing we attempt to constrain the luminosity at the epoch of 
emergence, using both analytic estimates and numerical simu- 
lations. 

4.2. Analytic Estimates for SN1997D 

In an analytic approach, we can use the average quantities 
of the helium-rich layer to estimate the radioactive heating and 
accretion luminosities, since this layer is expected to be both 
the source of material for late time fallback and the main con- 
tributor to 7-ray trapping. As a reference composition of the 
bound material we use the best fit model of Young et al. (1998 ) 
(see below) which is - by mass - H:0.10, He:0.45, 0:0.45. The 
electron scattering opacity at full ionization is n es = 0.22 and 
the atomic weight per electron is fi « 1 .27. 



4.2.1. Initial Density and Expansion Time Scale 

Note that for initially homologous expansion, the total ki- 
netic energy of the layer can be related to its total mass, M, and 
radius, R, according to 



E kb , = —MV" = —M\ ~ 



10 



10 



to 



(24) 



where V is the velocity at the outer radius and ?o =R/V is the 
corresponding expansion time scale. By expressing the radius 
in terms of the mass and initial average density, po we can write 



3/2 



M 



3/2 



(25) 



Equation ( p4| ) is especially convenient, since the average en- 
ergy per unit mass deposited in a roughly homogeneous layer is 
relatively easily determined in simulations (see Woosley 1988 
for SN1987A). The fraction of the total explosion energy car- 
ried as kinetic energy of the helium-rich layer is typically a 
few percent. Adopt ing the global quantities of the model of 



uratto et al. (1998|), we have in the case of of SN1997D a 

1.2 xlO 49 ' 



helium-rich layer with Mn e ~ 5 M, 
' ergs gm" 



'0, Eki„ He 

or Ekinfle/ M[ie ~ 1 .2 x 10 15 ergs gm" 1 . For this layer 



p f() ~ 9.4 x 10 9 gm cm 3 s 3 



ergs, 



(26) 



Note that the relatively low energy of the explosion is clearly 
reflected in th is quantity, wh ich for SN1987A is only 



10 9 gm cm" 3 s 3 (Chevalier, 1989|). 



The value of pqIq is sufficient for estimating 7-ray trap- 
ping (eq. [|T1||), but additional information is required for es- 
timating the accretion luminosity, which scales as (pQtffil^t^ 4 
(eq. This difference is an inevitable consequence of the 

fact that the bound material is not actually expanding homolo- 
gously, since the gravitational deceleration cannot be neglected. 
For a realistic supernova envelope there is no time instant when 
the entire envelope is in coherent homologous motion. As a 
first approximation, we assume that the entire envelope has a 
common expansion time, which turns out to be fo « 30 hrs for 
an initial radius of 2 x 10 13 and a hydrogen envelope of 18 M@. 
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We also have pHe,o « 7.45 x 10 6 gm cm 3 for these parameters. 
The total mass bound to a black hole of 3 Mq is then 

4-7T , 8-7T , 

M hnd = j/)ft,o«rt,D = —GM BH 4po. He = 0.136 M Q , (27) 
where we have used equation (p^). 

4.2.2. Accretion Time Scale 

Homologous expansion implies that the thermal energy is 
sufficiently smaller than the kinetic energy. In particular, in 
order to have E,h/Eki„ <^0.1 for the helium-rich layer presented 

above, the temperature in the layer must be To.He <> 8 x 10 5 °K. 
The corresponding gas and radiation pressure are of similar 
magnitude, leading to an estimate of the initial sound speed 
in this layer of c s .o(He) ss 10 7 cm s . The initial accretion 
time is thus t acc> o « 120 hours, significantly longer than the 
expansion time scale. In this analytic estimate, the accre- 
tion flow is approximately dust-like from its onset. As long 
as we assume that the total thermal energy is the order of 
~ 10% of the kinetic energy of the helium layer, which itself 
is well constrained, the temperature in this layer is approxi- 
mately fixed at the above value. The initial gas pressure scales 
as Po ~ Tpo, so the sound speed cannot decrease significantly 
below c Si o(He) w 10 7 cm s _1 (if the pressure is gas dominated 
Cj (He) oc Po/ po oc T), and the initial accretion time scale ( pj| ) 
cannot deviate sharply from t acc $ m 100 hrs, even if the combi- 
nation of po and fo is allowed to vary. 

4.2.3. Critical Time 

A key feature in this analytic estimate is that t cr u ~ 180 hrs 
(eq. ffijft) or ~ 6fo- As discussed in § ^, dust-like accretion 
should require a time of about to to settle into the late-time tem- 
poral decay. Hence it is likely that radiation forces will impose 
some modification to the accretion flow. We can expect that 
when the accretion flow finally resettles on a dust-like motion 
the total bound mass will be reduced, and the average density 
and expansion time in the bound region will be modified. The 
initial values of po and fo can thus provide an upper limit on the 
late time accretion rate and luminosity through equations ( p"6j ) 
and @. 

The effect of radiation pressure on the early-time accretion 
flow when the luminosity is close to the Eddington limit is non- 
linear, and therefore cannot be ascertained from our analytic 
treatment. However, by making two simplifying assumptions 
we can roughly estimate the impact of the Eddington-limited 
stage on the late-time accretion rate. First, we assume that the 
accretion flow adjusts so that the accretion luminosity is ex- 
actly the Eddington value. Effective gravity is then negligible 
everywhere. Thus, while the average density po(t) and expan- 
sion time to(t) maintain a constant po(O f o(f) 3 = Po f O' individu- 
ally they vary according to 



Po(t) = 



Rp + Vpt 
Rp 



The combination ppt^ 3 , which determines the magnitude of the 
dust-like accretion rate (eq. Jl6|]) evolves as 



fo(0 = k 



Ro + Vpt 
Ro 



(28) 



Po(f)^ /3 (0 = A)fo /3 ( 1 + f 



-1/3 



The second assumption is that the transition to a dust-like flow 
is instantaneous, and that at the time of transition, td ust , the val- 
ues p~o(t du st)wd %{t d ust) exactly yield M(t dusl ) = M crU through 
equation (jig). From equation (28) we find that tdust must then 
satisfy the relation: 



41.15 (JL^fJL 
\Q.5J V0.4 



1/2 / M, 



Mr. 



-1/3 



PO.He 



10" 4 gm cm -3 

-5/18 



5/6 



1 + 



tdust 

to 



tdust 

to 



( h).He \ 

Vlhr/ 



-25/18 



5/6 



1 . 



(30) 



(29) 



The modification due to radiation pressure is that tdust < 
t cr it, so that the late-time accretion luminosity is L(t) = 
LEdd(t/td us tT 25/l& instead of L(t) = L Edd (t /t C ri,T 25 ' 1S - We find 
that tdust ~ 129 hrs, so the late time accretion luminosity will be 
reduced bya factor of ~ 0.633 compared to that predicted by 
equation (p~8]) for the initial parameters. 

Note that in the limit tdust ^ k) equation ( [30| ) reduces to the 
form fj ((5/ oc potQ, so t dm , is dependent only on the initial spe- 
cific energy of the bound material. This is a natural conse- 
quence of the approximation used here, since ongoing homol- 
ogous expansion eventually loses memory of the initial size Rq 
(which determines to and po separately). For the values pre- 
sented above, tdust is larger than to by only a factor of a few so 
the above argument is not exact, but if fo is confined to a reason- 
ably small range, the late-time accretion luminosity can still be 
very well constrained by the specific energy of the helium-rich 
layer. Indeed, we find that t dml varies by less than 10% when fo 
is varied in the range 20-40 hrs if pofle^o is kept fixed. 



4.2.4. Radioactive Heating and an Estimate of Emergence 

With all the necessary values for estimating the late- 
time accretion rate and 7-ray opacities, we now proceed 
to examine the competition between the accretion and ra- 
dioactive luminosities. The exact time and luminosity at 
emergence depend on the abundances of 57 Co and 44 Ti in 
the envelope. These two isotopes are mostly produced in 
deeper layers than those where 56 Co is synthesized (espe- 
cially 44 Ti , which is produced through Nuclear-Statistical- 
Equilibrium (NSE) in the innermost layers of the shocked 
envelope material, Timmes et al. 1996). We examine their 
effects with two extreme assumptions - either that they are 
completely absent, or that their abundances scale with that of 
56 Co . In the latter case, using M 5 7 Co /Ms6 Co (SN1987A) = 
0.0243, M 4 4 Ti /Ms6 Co (SN1987A) = 1.33 x 10" 3 we obtain 
M 57Co (SN1997D) £ 5 x 10" 5 M 0) M« Ti (SN1997D) < 2.7 x 
10~ 6 Mq. Assuming that the opacity to 7-rays is indepen- 
dent of the accretion history (as is reasonable, since the initially 
bound region is a small fraction of the total helium-rich layer), 
we can treat the radioactive heating history as fixed once the 
initial combination po,Heto is given (eq. JlO|]). 

The resulting luminosities due to accretion and radioactive 
heating for our reference parameters of fo = 30 hrs and po.He = 
7.45 x 10~ 6 gm cm~ 3 are shown in Fig. |l|. If the abundances of 
57 Co and 44 Ti are negligible, the total luminosity due to ra- 
dioactive heating is equal to the heating rate through decays 
of 56 Co . In this case, the accretion luminosity does emerge 
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Table 1 

Estimates of black hole emergence in SN1997D based on the analytic model (see text). 
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a Co and Ti abundances scaled according to SN1987A 
b No 57 Co and 44 Ti (only 56 Co ) 

' In this model the accretion luminosity does not reach 50% of the total luminosity (although it does constitute more 
than 40%) 



with L acc {tBH) = L56 Co (?bh) ~ 3.6 x 10 35 ergs s"'at a time of 
t BH ps 1216 days (April 2000), after which the light curve will 
settle on a power law decline rather quickly. If 57 Co and 44 Ti are 
present with abundances rescaled from SN1987A their impact 
on the luminosity at ~ 1000 days is not negligible, and accre- 
tion luminosity must be compared to the total radioactive heat- 
ing (labeled L, ot - rac i in Fig. nh. In this case emergence is de- 
layed to Ibh ~ 1466 days (December 2000), with L acc (tBH) = 
Ltot-md(tBH) & 2.9 x 10 35 ergs s" 1 . Note that in this case the con- 
tribution of 44 Ti is about 1.4 x 10 35 ergs s" 1 , which decreases 
only very slowly due to increasing 7-ray transparency. Cor- 
respondingly, the total luminosity (labeled L, or ) does not fol- 
low a perfect power law until 44 Ti heating becomes negligi- 
ble, although the accretion luminosity will remain the dominant 
source of luminosity for several thousands of days after tBH- 

The results of similar calculations for the emergence of a 
3 Mq black hole in SN1997D for different values of the ini- 
tial parameters are presented in table [j] We vary the energy of 
the 5 Mq helium-rich layer as 1, 3 (the nominal case) and 5 per- 
cent of the total explosion energy of 4 x 10 50 ergs, while vary- 
ing to in the range 20-40 hrs. We list the derived values of f m , 
and tdust for eacn model and the resulting luminosity and time 
of emergence. Due to uncertainty regarding the abundances of 
56 Co and 44 Ti , we examine the two extreme cases where their 
abundances scale according to SN1987A, and when they are 
completely absent. The case of the helium-rich layer having 
5% of the total energy and with the significant abundance of 
57 Co and 44 Ti marks the border-line of achieving emergence 
after a few years (instead of in thousands of years) as the the 
combination of 57 Co and 44 Ti heating is slightly larger than the 
accretion luminosity (by about ~ 20%) for a very long time. 

Our results indicate that if a black hole was formed, it is 
likely to emerge about three years after the explosion. Taking 
into account also for the uncertainties in composition and black 
hole mass, we estimate that plausible values of the luminosity at 
emergence seems to lie in the range 0.3-3 x 10 36 ergs s" 1 , and 
the expected time of emergence is ~ 1100- 1500 days after 



the explosion, with higher luminosities corresponding to ear- 
lier times. We emphasize that the dominant quantity that deter- 
mines the value of the luminosity at emergence is the kinetic en- 
ergy of the helium-rich layer, while the time of emergence and 
the character of the late-time light curve are also sensitive to 
the possible presence of 57 Co and 44 Ti in non-negligible abun- 
dances. 
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FIG. 1. — Analytic estimates for the total bolometric luminosity of 
SN1997D, L tot , and the partial contributions of accretion onto the black hole 
(L acc ) and of radioactive heating (Lt„t-rad = Ls6 Co +Lsi Co +Ln Ti ) during the 
period 800- 1600 days from the explosion (Jan. 1999-Mar. 2001). Arrows 
mark the time of emergence of the accretion luminosity in the case where ( 1 ) 
large arrow: s7 Co and 44 Ti are assumed to be present in the envelope with 
abundances scaled down from SN1987A (see text) and (2) small arrow: when 
only 56 Co is present. 
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4.3. Numerical Results for SN1997D 

While an approximate analytic treatment and scaling behav- 
ior can provide useful estimates of the time and luminosity at 
emergence, it also underscores the need to perform numerical 
investigation which can track in detail the early time accretion 
history and the late-time evolution of the light curve. The nu- 
merical simulation is also required to quantitatively assess the 
role of non-linear radiation-hydrodynamics processes and the 
effects of variable chemical composition, realistic opacities and 
radioactive heating. In order to explore in detail the emergence 
of a black hole in SN1997D, a full radiation-hydrodynamic sim- 
ulation of fallback of material from the supernova envelope up 
to the several years after the explosion has been performed, us- 
ing the numerical code described in § 3. The initial conditions 
for the ejecta were based on the best-fitted model of Y oung et 
al. (1998) for the SN1997D explosion. This model of a 26 M 
progenitor was fitted on the base of the early time light curve 
(up to ~ 100 days after the explosion), which is practically in- 
dependent of the inner part of the helium-rich layer. Hence, 
while a simulation based on this model serves as an essential 
quantitative estimate, the unavoidable uncertainties regarding 
the initial velocity, density and temperature profiles of this in- 
ner layer obviously limit the applicability of the simulation as a 
"best fit" also for the accretion history and luminosity. 

4.3.1. Initial Profile 

The initial profile describes the supernova envelope at the 
time of break-out of the shock at the surface, corresponding to 
about 13 hours after core bounce. The expanding envelope is 
composed of two main components: an inner helium-rich layer 
and an outer hydrogen-rich envelope. The properties of the pro- 
file are presented in Figs. || and 3] and the region of initially 
bound material is shown in more detail in Fig. By assuming 
that the abundance of radioactive isotopes scales with the oxy- 
gen mass fraction, we can simply estimate the mass fraction of 
56 Ni in every mass shell by requiring that the total 56 Ni mass 
is 2 x 10~ 3 Mq (equal to the total 56 Co mass inferred by obser- 
vations). In the simulation we assume that the abundances of 
57 Co and 44 Ti scale with the abundance of 56 Ni , so their total 
masses are set according to the aforementioned quantities. 

In the context of late-time accretion, the most dominant fea- 
ture of the initial profile is the non-uniformity of the helium-rich 
layer. A significant fraction of the mass is concentrated in an 
over-dense region bordering with the hydrogen-rich layer. This 
region also holds the bulk of the kinetic energy of the helium- 
rich layer, 1.7 x 10 15 ergs gm" 1 , somewhat larger than the aver- 
age value of 1 .2 x 10 15 ergs gm" 1 mentioned above. The bound 
region on the other hand, is relatively under-dense and slow. 
The total initially bound mass is ~ 0.22 Mq. 

The expansion time scale throughout most of the initially 
bound material is fo ~ 50-60 hrs. This region is radiation- 
pressure dominated, and the initial sound speed is roughly 
c s .o(He) « 1.3 x 10 7 cm s" 1 . The corresponding accretion time 
scale, t acc x) ~ 50 hrs is therefore similar to the initial expansion 
time scale. Unlike the analytic cases discussed above, in this 
model pressure forces are initially important. The long initial 
expansion time scale more than compensates for the decreased 
density in terms of determining the accretion rate, and we have 
f„„ ~ 300 hrs. This model is also expected to reach very high 
accretion rates at early times, which will result in some mod- 
eration of the accretion flow when M cr - U m 3. 13 Mq yr is ap- 
proached. 
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4.3.2. Computational Aspects 

Tracking the evolution of the envelope over several years (the 
expected time of emergence of the black hole) is not feasible 
numerically. Integrating the model until emergence would re- 
quire a CPU time greater than one month, even when adjusting 
the position of the inner boundary and using the MTP proce- 
dure (see § ||). As a result, we devised a rescaling scheme where 
the original model is rescaled to a more rapidly evolving one, 
while maintaining simple relations between the quantities com- 
puted for the rescaled model and the original one. This rescal- 
ing scheme, presented in detail in the appendix, has allowed for 
an additional acceleration of a factor of ~ 5 in computational 
efficiency, thus reducing the total required computational time 
to an acceptable limit. 

We must note that several particulars of the initial model used 
here for SN1997D imposed some specific complications in em- 
ploying the rescaling scheme. The fact that initially t acc < to and 
that the flow is expected to be moderated at early times when 
a critical accretion rate is approached implies that the rescaling 
scheme described in the appendix cannot be applied at the ear- 
liest stages of evolution. Furthermore, at early times the energy 
generation rate of the radioactive elements, and especially 56 Ni , 
is at its largest. The enhancement required for these rates in the 
course of rescaling (see in the appendix) further compounds the 
applicability of the rescaling scheme at very early times. 

Therefore our numerical approach has been to perform the 
simulation of the model in its original scale until three condi- 
tions are met: 1) the time dependent marginally bound radius 
is significantly smaller (a factor of 5) than the accretion radius; 
2) the accretion rate through the inner boundary is smaller than 
Q.\M cr it and 3) that total energy to be emitted in the decays of 
the remaining 56 Ni is significantly smaller than the internal en- 
ergy at that time. We find that these conditions are all satisfied 
at the physical time of about t =25 days, which is when the sim- 



ulation was stopped, rescaled by a factor of 5, and restarted. We 
performed the simulation of the rescaled model with a single 
time-step integration until the recombination front has settled 
to a roughly constant Lagrangian position, which occurs after 
about t = 150 days. These first two stages required about 20 
CPU hours each. From this point on, the simulation is contin- 
ued with the Multi Time-step Procedure, where the integration 
domain was divided into four sub-grids. Performing the simula- 
tion up to a physical time of t = 1500 days required an additional 
120 CPU hours (recall that at later times the time step drops as 
the inner boundary must be moved inward in order to maintain 
it in LTE - see § ||). All computations were performed on a SGI 
500MHz EV6 machine with a Compaq XP1000 processor. 

4.3.3. Numerical Results 

The calculated light curve for SN1997D is shown in Fig. |[ 
Also shown are the observed data exten ding up to ~ 440 days 
after the explosion ( Benetti et al., 200C| ), with which the agree- 
ment of the simulated light curve is very good. For comparison, 
we also show the light curve calculated for an identical model 
but with the abundance of radioactive elements reduced by a 
factor of 300, and also the late-time light curve estimated ana- 
lytically as in Fig. [j] above. In the simulation we cannot single- 
out the accretion luminosity directly, by we can evaluate its con- 
tribution to the total luminosity by subtracting the contributions 
of radioactive heating. The simulation does record the total rate 
of energy deposition in the envelope by the decays of each of 
the isotopes, Qxif), and since we can estimate that the deposited 
heat is emitted as a bolometric luminosity very rapidly after 
thermalization, i.e., Lx(t) ~ Qx(t), we can attempt to identify 
the accretion luminosity as L acc = L tot - (Q^ Co + 2 57 co + 2 44 ri)- 
These calculated contributions for the model are shown in 
Fig-! 
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The transition from a Co -dominated tail in the luminos- 
ity is evident at about 1000 days. From Fig. || we can de- 
termine that the accretion luminosity does become the dom- 
inant source of the light curve, and the time of emergence 
(when L acc (t BH ) = \L tot (t BH )) is t BH rss 1050 days, when L acc « 
3.2 x 10 35 . This is roughly i of the accretion luminosity (at 
the same physical time) estimated using the initial average val- 
ues of po t He = 5 x 10~ 6 gm cm~ 3 and to = 50 hrs in equation (p"8]). 
This difference is the result of moderation by radiation forces 
as the accretion rate exceeds M cnf , and in part also by the fact 
that initially the envelope's own pressure is not negligible. An 
accretion rate of M > M cr - U is actually reached while the accre- 
tion flow is still building, and the radiation pressure it induces 
causes a significant part of the initially bound material to be- 
come unbound. Indeed, the the total accreted mass (extrapo- 
lated to t —> oo) found in the simulation is decreased to only 
~ 0.067 M©, less than a third of the originally bound material. 

The evolution of the accretion rate is plotted in Fig. [7J, which 
shows the accretion rate (calculated near the inner boundary) as 
a function of time. The maximum accretion rate is reached af- 
ter about 18 hrs, when its value is ~ 1.5M cr!f . At this stage the 
accretion flow is not quasi-stationary (see below), and signifi- 
cant moderation is imposed. We comment that the details of the 
accretion rate at these very early times are dependent on the as- 
sumed initial conditions. With initial conditions set during col- 
lapse, MacFadyen, Woosley, & Heger (1999) find that the earli- 
est accretion reaches a maximum rate as early as a few hundred 
seconds after bounce. In our numerical simulation the initial 
conditions are set after the passage of the shock through the en- 
velope and hence we cannot reproduce these earliest stages of 
the accretion history. More significant in the context of black 
hole emergence, however, is the character of the flow as it set- 
tles on the late-time, dust-accretion track, where there is good 
agree ment between our results and those of M acFadyen et al. 
(1999). 
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FIG. 6. — The total bolometric luminosity for the nominal model for 
SN1997D and the estimated partial contributions of accretion onto the black 
hole and of radioactive heating during over the period 700- 1500 days from 
the explosion (Nov. 1998-Dec. 2000). The arrow marks the time of emergence 
of the accretion luminosity (L acc = jL,ot)- 
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to the rescaled one (see text). 



It is apparent that the accretion does eventually settle on a 
dust-like flow M~f" 5 / 3 , as expected, at ~ 10 days after the 
explosion. The quantitative effect of the early history of the ac- 
cretion flow is clearly apparent, since the accretion rate is lower 
than the value predicted by equation ( p"o| ) forM(10 days), based 
on the initial properties of the inner part of the bound region. 
By extrapolating the dust-like accretion rate back in time, we 
find that M = M cr u at ~ 100 hrs, instead of the original estimate 
of t cr it ~ 300 hrs. 

We mark the time of transition between the original model 
and the rescaled one by a circle in Fig. ^. After a short tran- 
sient, the accretion flow resettles on a dust-like accretion flow, 
with a very good fit (accurate to about 10%) to the initial un- 
sealed model. This agreement suggests that the rescaled model 
is indeed capable of reproducing the accretion flow and lumi- 
nosity. 

In Fig. H we show the behavior of the computed accretion 
luminosity as a function of the computed accretion rate M, 
and compare it to the analytic estimate of the Blondin formula 
(eq. [|l7]]). The excellent fit of the simulation to the analytic ra- 
tio demonstrates that the radiation-hydrodynamic evolution in 
the accreting region does indeed follow a sequence of quasista- 
tionary states, so the analytic estimates are also justified in the 
case of a variable chemical composition. The reduced accretion 
rate is clearly the cause for the accretion luminosity being only 
about one fifth of its value (had the entire initially bound mass 
remained available for accretion). 

It is noteworthy that the black hole emerges somewhat earlier 
than in the case of Fig. [j], in spite of the lower accretion lumi- 
nosity. This is due to a significant complementary effect: the 
structure of the helium rich layer reduces the effective optical 
depth to the high energy photons emitted in the radioactive de- 
cays. First, the initial column depth of the layer is only about | 
of that of a layer with identical mass and size but a constant den- 
sity. Second, the outer part of this layer has a larger specific ki- 
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netic ene rgy and therefore a smaller pot^ than the average value 
used in § 4.2. This combination causes the 7— ray column depth 
of the helium-rich layer to be only 1/5-1/4 of that at an iden- 
tical time for an analytic estimate, and the trapping efficiency 
is decreased accordingly. Correspondingly, the contribution of 
radioactive heating to the bolometric luminosity is reduced to 
what would have been expected if the entire helium-rich layer 
were characterized with po = 5x 10~ 6 gm cm" 3 and to = 50 hrs. 

The model with reduced radioactive isotopes was investi- 
gated in order to assess the quantitative effect of radioactive 
heating on the dynamics of the flow. This effects can be inferred 
indirectly from Fig. Bl and directly by comparing the accretion 
history in both models. We find that the accretion histories are 
practically identical, where the late-time accretion rate in the 
model with reduced radioactive heating is larger by only ~ 5% 
than in the nominal model, implying that even in the nominal 
model the energy deposited by heating hardly affects the region 
of bound material. The calculated accretion luminosities also 
agree to within a few percent, although it is the model with ra- 
dioactive isotopes that suggests a higher luminosity. While this 
effect may be partially a numerical artifact, it may also arise 
from the slightly larger photo-spheric radius found in the origi- 
nal model, presumably caused by radioactive heating. We con- 
clude that when the abundance of radioactive elements in the 
ejected envelope is as low as it is in the case of SN1997D, it 
is unlikely to impose a significant effect on the expected accre- 
tion rate or luminosity. By contrast, even the small amount of 
radioactive isotopes present in the SN1997D envelope is suffi- 
cient to affect the decline of the recombination peak (see Fig. |]). 
This effect was demon strated clearly in the investigations of the 
SN1 987A light curve (|Woosley, 19881; Shigeyama & Nomoto, 
1990), and is mainly due to release of heat deposited in the in- 
ner part of the envelope by 56 Ni decays. 



35.8 



♦ ♦ ♦ 



35.6 



0J 



SO 
O 



♦ ♦ 



35.2 



35.0 



-4.0 



-3.8 



-3.6 



-3.4 



M (M a /yr) 



FIG. 8. — The computed accretion luminosity (Fig. ^) versus M (diamonds) 
for the fiducialjiiodel of SN1997D. Also shown is the analytic estimate based 
on equation (JlVt) for the chemical composition of the bound region of the 
helium-rich layer (mass fractions H:0. 1, He:0.45, 0:0.45). 

The negligible effect of radioactive heating on the hydrody- 
namic evolution allows for a simple extrapolation regarding the 
expected light curve if the abundance of the relevant isotopes is 



varied. If the abundances of 57 Co and **Ti are very low, emer- 
gence can be gauged by comparing the accretion luminosity 
and the heating by 56 Co decays, suggesting that emergence oc- 
curs at about 970 days. The luminosity at emergence would be 
slightly higher, with L acc m 3.9 x 10 35 . In this case, the contri- 
bution of 56 Co heating declines rapidly and the late time light 
curve should indeed fall off as a power-law in time. Note that a 
perfect power law decline is found, as expected, for the model 
with reduced abundances of radioactive isotopes. 

To conclude our investigation of the SN1997D light curve, 
we examine the history of the luminosity and the accretion flow 
throughout the envelope. The profiles of the (comoving frame) 
luminosity, L(r), and mass flux, |M|(r), at various selected times 
during this evolution are plotted in Figs, ^a and^i respectively. 
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The luminosity profiles demonstrate the evolution of the light 
curve through its various stages (see Zampieri et al. 1998a for 
details). During the stage dominated by radioactive heating the 
photo-sphere lies well above the accreting region. The presence 
of a separate source of luminosity due to compressional heat- 
ing in the course of accretion is evident in the inner part of the 
flow, but only when its magnitude becomes comparable to the 
heating produced by radioactive decays does it start to dom- 
inate the luminosity output. The drop in luminosity near the 
inner boundary arises because the opacity of the partially ion- 
ized helium-rich material becomes larger (much greater than 
the electron scattering opacity). Closer to the black hole the 
increased temperature increases the degree of ionization of the 
helium-rich material, and the comoving luminosity rises again 
to a local maxima. Note that the mass flux remains self sim- 
ilar through out the evolution, as the accreting region grows 
in size but the expansion of the envelope reduces the magni- 
tude of the accretion rate. The sharper, inner spike in the mass 
flow rate is due to the original spike in density of the model at 
the interface between the helium-rich and hydrogen-rich layers 
(the smoother, outer peak simply corresponds to a maximum 
in the combination pvr 2 ). In the inner region the accretion rate 
is nearly independent of radius, which also reflects that the ac- 
cretion flow does indeed follow a sequence of quasi-stationary 
states. 



5. BLACK HOLE EMERGENCE IN OTHER SUPERNOVAE 

While SN1997D appears to be the only existing candidate 
for observing black hole emergence, theory suggests that the 
explosions of the most massive stars in the core-collapse su- 
pernovae range would be significantly more favorable for an 
unequivocal determination of the presence of a black hole. The 
key feature is, of course, the final yield of heavy elements in 
the envelope following early fallback. These heavy elements 
include the radioactive isotopes which mask the accretion lu- 
minosity for extended times, so that the lower the abundance of 
radioactive elements, the shorter the time until emergence and 
the larger the luminosity at that time. SN1997D may be suffi- 
ciently poor in heavy elements so that the accretion luminosity 
can dominate the heating by 44 Ti decays, but 56 Co decays are 
by far the dominant source of power until ~ 1000 days after the 
explosion. 

In a survey o f the supernovae of massive stars, W oosley 
& Weaver (1995) concluded that explosions of stars of initial 
masses of 30 - 40 Mq are likely to include significant early 
fallback, thus advecting the newly synthesized heavy elements 
onto the remnant, a likely black hole. For "nominal" explo- 
sion energies of ~ 1.2x 10 51 ergs, they found that explosions 
of stars with masses of 30-40M Q should easily lead to rem- 
nant masses in the range of 3- 10M Q , while ejecting envelopes 
practically free of heavy radioactive elements. While they did 
not consider explosion energies as low as the 4 x 10 50 ergs in- 
ferred for SN1997D, it is natural to assume that lower energy 
explosions will give rise to enhanced earlier fallback and thus 
higher remnant masses 8 . 

In order to gauge the possible signature of black hole forma- 
tion in a supernova of a very massive progenitor, we computed 
the light curve of such a supernova using the numerical code 
described in § [| As a reference we used a 35 Mq progenitor 



based on model S35A of Woosley & Weaver (1995 ). It has a 
pre-explosion outer radius of 8 x 10 IJ cm, and is composed of 
a 14 Mq helium-rich core and a 21 Mq hydrogen-rich enve- 
lope. The post explosion ejecta profile was constructed so that 
the total ejected masses are 11.5, 12 and 4 Mq of hydrogen, 
helium and oxygen, respectively. Most notably, the envelope 
is assumed to be free of radioactive isotopes. The total kinetic 
energy of the ejecta is set at 1.2 x 10 51 ergs, leaving behind a 
7.5 Mq black hole. Assuming that the helium-rich layer car- 
ries 1% of the total energy and that the entire envelope has 
a common expansion time scale, our model has to ~ 55 hrs. 
The initial profile of the bound material thus has po^e = 4.93 x 
10~ 6 gmcm~ 3 and therefore po?o = 3.83 x 10 10 gmcm~ 3 s -1 . This 
expansion time scale is significantly shorter than the accretion 
time scale, which for this model is roughly 10 days, so that 
the flow is dust-like at its onset. The critical time is, however, 
~ 280 hrs, and once again we expect the accretion can begin 
as dust-like and become then moderated by radiation pressure 
when the Eddington limit is approached. Through the analytic 
estimates in § \.2 we find that tdust « 110 hrs for this model. 



The calculated light curve for S35A is shown in Fig. 10 We 
note that no rescaling was required for these models due to the 
relatively large time step and short evolutionary time required 
until emergence; the calculation required about 55 CPU hours. 
The modest explosion energy combined with the large radius of 
the progenitor lead to a relatively small initial average temper- 
ature of about 2.5 x 10 s °K, and recombination occurs around a 
time of t rec w 35 days. The light curve during this peak is quite 
similar to the one observed for SN1994W, a reasonable can- 
didate for an explosion of a very massive star with very little 
radioactive elements in the ejecta (see below). 
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8 Woosley & Weaver (1995) did examine explosion energies of > 2x 10 51 ergs 



Fin 1 n — Thf light mrvp inHiiHincr an accretion tail for the S35A model 
of Woosley & Weaver (1995). Solid line: Mbh = 7-5 Mq, dashed line: 
Mg H = 15 Mq. 

In the absence of radioactive heating the recombination peak 
fades rapidly, and the emergence of the black hole occurs at 
f « 55 days with a total luminosity at emergence of about 

In general, they found that if explosion energy is large enough to eject a non- 



negligible amount of radioactive isotopes in stars of 30—40 Mq, the remnant mass will be < 2.5 Mq, and hence, possibly, a neutron star 
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10 37 ergs s -1 . The rapid decline of the recombination luminos- 
ity causes the accretion luminosity to be nearly the sole contrib- 
utor to the light curve as early as at t «60 days, with a luminos- 
ity of L(t = 60 days) w 7.2 x 10 36 ergs s" 1 . This value is smaller 
than the analytic result of L(t) = L£j I ;(60days/f £ ; i(Jf )" 25 / 18 by a 
factor of a few, although the calculated accretion rate does sat- 
isfy M(t) ~ M C rit(t ftdustT 5 ^- This result is due to the fact that 
at this stage bound-bound and bound-free opacities in the (rel- 
atively narrow) partially ionized region in the accreting mate- 
rial are high (much larger than free electron scattering opacity). 
This imposes a significant contribution to the total optical depth 
of the accretion flow, increasing its "effective" average opacity. 
Note that for a given accretion rate, the accretion luminosity is 

1/2 



expected to follow ' dependence (using eq. [Ill] 

and expressing Lgdd and Msdd explicitly). This is a transient 
effect: as time elapses, expansion of the envelope gradually 
causes the fully ionized, inner part of the accreting region to 
dominate the optical depth of the flow, so that the average "ef- 
fective" opacity is rapidly declining towards the electron scat- 
tering value. Consequently, at this stage the luminosity falls 
off somewhat more slowly than a f~ 25 / 18 power-law, eventually 
settling on such a time dependence. 

Also shown in Fig. |o|is the same model with the black hole 
mass arbitrarily increased to Mbh =15 M Q , the light curve evo- 
lution is identical until emergence, which occurs at a slightly 
higher luminosity and therefore at a slightly earlier time (~ 3 x 
10 37 ergs s~'atf « 50 days). The light curve again evolves there- 
after due to accretion L(60days) « 1.33 x 10 37 ergs s" 1 , again 
falling off somewhat more moderately than a f~ 25 / 18 power law. 

We emphasize that the accretion rate and luminosity in both 
models are significantly lower than that predicted by equa- 
tions (|l6[ [l^) for the initial parameters of the helium-rich layer. 
Even though the accretion can begin as dust-like, the modera- 
tion of the accretion flow by radiation pressure when the Ed- 
dington limit is approached is evident in the absolute magni- 
tude of the accretion luminosity-tail. The model with the larger 
black hole mass has a higher critical accretion rate, so even 
though it also induces a smaller tdust, the ra ti° of the accretion 
luminosity is increased by slightly more than just a factor of 

2 2/3 

arising from the ratio of the black hole masses. 
We conclude, that Type II supernovae with no radioactive 
isotopes provide the most favorable cases for observing the 
emergence of a black hole. It is noteworthy that a similar ex- 
plosion of the bare helium core in Type Voile supernovae would 
be far less favorable. Not only will all the explosion energy 
be deposited in the helium layer (instead of only a few per- 
cent), there will also be no inverse shock that slows the in- 
ner part of the helium layer and enhances fallback (which is 
essential for decreasing the abundance of ejected radioactive 
elements). The very small time step required for simulating 
such a model prevents us from examining it numerically, but 
we can place some limits with the analytic derivations pre- 
sented in § \.2. For example, consider a 14 M Q helium star 

10" cm exploding with an en- 



with an initial radius of Rq 

ergy 1.2x 10 51 ergs, producing, again, a 7.5 M black hole. We 
estimate that the post-explosion envelope will have po.Hefo ~ 
1.85 x 10 7 ergs gm" 1 , po ~ 3 gm cm" 3 and to w 0.05 hrs. Even 
assuming that the flow does begin as dust-like, we find that 
tdust ~ 4 hrs (t cr i t k 10 hrs) so that the luminosity at t ~ 60 days 
will be only ~ 5.3 x 10 35 ergs s" 1 - a few percent of the values 
found above. While recombination is likely to occur earlier in 
such supernova, the actual accretion rate in such an explosion 



will be strongly modified due to high initial temperatures (and 
pressure), confirming that even radioactive-free Type Ib/Ic su- 
pernovae are unlikely to be good candidates for observing the 
emergence of black holes. 

6. CONCLUSIONS, DISCUSSION AND OBSERVATIONAL PROSPECTS 

In this work we have examined the emergence of luminosity 
due to spherical accretion onto black holes formed in realistic 
supernovae. The study is based on a combination of analytic 
estimates and numerical simulations. Our main focus has been 
on SN1997D, where the low explosion energy and very low in- 
ferred abundance of 56 Co in its ejecta provide a possible (and 
currently, the only) candidate where the emergence of a black 
hole in a supernova may actually be observable. 

We confirmed the mai n features of the accret ion history and 
luminosity as derived by Zampieri et al. (1998c). At late times, 
the accretion flow settles on a dust-li ke motion, and the accre- 
tion rate falls off as a ?" 5 / 3 in time (Z!olpi et al., 1996). Ac- 



cretion proceeds in a sequence of quasistationary states, so the 
accretion luminosity achieves the ex pected magnitud e of a hy- 
percritical, spherical accretion flow (Blondin, 1986), account- 
ing for the secular decay in time of the accretion rate. The 
resulting time dependence of the late-time accretion luminos- 
ity is L(t) oc f" 25 / 18 . This relation marks a fundamental feature 
regarding emergence: since the main competitor to accretion in 
powering the late time light curve is radioactive heating which 
decays exponentially with time, emergence of the accretion lu- 
minosity is inevitable if spherical accretion does indeed persist. 
While in typical Type II supernovae the abundance of radioac- 
tive isotopes is expected to be large enough to prevent a practi- 
cal observation of the emergence, more rare explosions where 
this abundance is significantly reduced and fallback is enhanced 
provide important exceptions. 

The case of SN1997D o ffers a unique opportun ity to exam- 
ine black hole emergence (Zampieri et al., 1998b), in view of 



its inferred explosion energy of ~ 4 x 10 3U ergs s" 1 and very 
low abundance of 56 Co , i.e., only 2 x 10~ 3 M P) (T uratto et al., 
1998; penetti et al., 2000| ). We have have examined the ex- 
pected emergence of a black hole in SN1997D analytically and 
numerically. In the analytic study we noted that the early-time 
accretion is almost certain to generate an Eddington-rate lumi- 
nosity, which will modify the accretion history and decrease the 
late-time accretion rate. In the numerical investigation we have 
included the effects of a variable envelope composition. In par- 
ticular, we incorporated heating due to radioactive decays, in- 
cluding a finite (time-dependent) 7-ray optical depth. The wide 
diversity of time scales involved compelled us to use a rescaling 
scheme in the numerical simulations (see in the appendix). 

Our various analyses offer a consistent assessment of the 
emergence of the blac k hole in SN1997D. We confirm the 
preliminary estimate of Zampieri et al. (1998b) that a 3 M Q 
should emerge in SN1997D about 1000 days after the explo- 
sion, in contrast to hundreds of years in the case of "stan- 
dard" Type II supernovae (like SN1987A) where the abun- 
dances of radioactive elements are significantly higher. We 
show that the dominant parameter in determining the lumi- 
nosity at emergence is the amount of kinetic energy carried 
in the post-explosion helium-rich layer (which is the source 
of late-time accretion). For realistic values of this energy, 
we estimate that the total luminosity at emergence will lie 
in the range 5 x 10 35 <J L to , <; 3 x 10 36 ergs s -1 , emerging at 
1000- 1500 days after the explosion. In our numerical study 
we have examined the best fit model of Turatto et al. (1998) 
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for SN1997D. The simulations suggest that, for this particu- 
lar model, the black hole emerges at about 1050 days after the 
explosion, with L ro , w 6.5 x 10 ergs s -1 . The results of the 
simulation, which is based on a consistent, fully-relativistic, 
radiation-hydrodynamics code, provides significant support to 
the analytic approach. Most notably, it demonstrates that, after 
passing through an early radiation-pressure-limited stage which 
reduces the accretion rate, the accretion flow does settle on the 
dust-like solution and that the influence of the radioactive de- 
cays on the hydrodynamic evolution is negligible for the low 
abundances inferred in SN1997D. 

The time of emergence and the corresponding luminosity de- 
pend on the abundances of 57 Co and 44 Ti in the envelope, which 
are unknown a-priori. However, as long as the accretion lumi- 
nosity is significantly larger than the threshold set by positron 
emission in 44 Ti decays, both emergence time and luminosity 
will not vary substantially. At the time of emergence 7-ray es- 
cape from the envelope is non-negligible and increases rather 
rapidly with time, as the envelope's optical depth drops. Corre- 
spondingly, the bolometric luminosity due to radioactive heat- 
ing decreases even faster than exponentially with time. 

Although our results appear to be quite robust within the con- 
text of the present study, some uncertainties arise due to effects 
not considered here. Most critical is the assumption of perfectly 
spherical accretion: if the bound material can maintain suffi- 
cient angular momentum, a centrifugal bar rier would drive the 
flow into a disk structure ( |Chevalier, 1996 ). The radiative effi- 
ciency of hypercritical disk accretion is uncertain, complicating 
any a-priori estimate regarding the character of the late-time ac- 
cretion luminosity. The residual angular momentum is likely to 
be dependent on the details of the progenitor and the explosion 
(we note that for SN1987A the obse rved neutrino bu rst seemed 
consistent with nonrotating models ( Burrows, 1988 ), indicating 
that angular momentum of matter close to the collapsed core did 
not play a significant role). In principle, late time accretion oc- 
curs many viscous time scales after the explosion, and much of 
the initial angular momentum could be lost, but the minimum 
specific angular momentum required to force disk accretion is 
quite small (about (GMbhRisco c ) 1 ' 2 > where Risco = 6GM/c 2 
is the innermost-stable circular orbit near a black hole). Fur- 
thermore, even if at late times angular momentum is negligible, 
some residual effect may arise if the character of the early-time 
accretion has been sufficiently altered (Mineshige et al. 1997 
found through simulations that sufficient angular momentum 
strongly affects the magnitude and time dependence of the ac- 
cretion during the first hours after collapse). We defer further 
investigation of fallback including angular momentum to future 
work. 

Additional effects could arise due to convection, mixing and 
clumping of envelope material (and radioactive isotopes in par- 
ticular), and dust extinction, although we expect that at a time 
of ;> 3 years after the explosion, any dust formed will have be- 
come transparent to the underlying luminosity. This may not be 
true if photoionization by the central accreting black hole (not 
included in our calculation) is important. 

Our results illuminate the uniqueness of SN1997D. It is the 
low abundance of radioactive elements, combined with the 
low explosion energy in general, that allows the black hole to 
emerge only a few years after the explosion, at which time the 
absolute magnitude of the accretion luminosity is still quite 
large (recall that for "standard" explosions we expect emer- 
gence luminosities of 10 32 - 10 33 ergs s" 1 ). Since gas and ra- 
diation are in LTE in the inner accreting region and the pho- 



tosphere is well localized in radius, the emerging spectrum is 
expected to be very similar to a blackbody at the photospheric 
temperature (~ 7500 K) with the superposition of heavy line 
elements. We therefore expect that emission will be mainly in 
the optical band. At a distance of ~ 14 Mpc, the correspond- 
ing apparent magnitude for the range of emergence luminosities 
predicted for SN1997D is m v w 28 - 30 and m R ~ m 7 « 27-29. 
The higher luminosity (lower magnitude) end of this range co- 
incides with the detection threshold of the HST STIS camera, 
indeed making emergence in SN1997D potentially observable. 
We note that the capabilities of HST are needed both because 
of the faintness of the source at this time, and the resolution, 
required to distinguish the supernova over the background of 
the parent galaxy NGC 1536. A recent estimate (M. Turatto 
1999, private communication) suggests that the required expo- 
sure time for the HST STIS to resolve the apparent magnitude 
expected at emergence in SN1997D at a signal to noise (S/N) 
ratio of 10 would be ~ 25000 seconds. The accretion lumi- 
nosity should be distinguishable from radioactive heating (or 
from circumstellar interaction, which has not been observed in 
SN1997D so far, Benetti et al. 2000) due to its unique power 
law dependence. This would be especially true if the abun- 
dances of 57 Co and 44 Ti are negligible, since after emergence 
the luminosity would be due only to accretion. But, even if 
there is a finite contribution from these two isotopes to the light 
curve, the presence of a power-law source may still be identi- 
fiable if several measurements with reasonable accuracy could 
be made. The novel consequences of a positive observation 
are self-evident - the first direct observational evidence that a 
black hole can be formed by an otherwise "successful" super- 
nova. Although our main focus has been SN1997D, other the- 
oretical Type II supernova models offer more favorable condi- 
tions for an actual observation of bl ack hole emergence. The 
survey of Woosley & Weaver (1995) suggests that progenitors 
with masses in range 30-40 M Q are quite likely to leave behind 
a remnant of mass 3 - 10 M and a practically radioactive-free 
envelope, as all the iron group isotopes are engulfed by the early 
fallback. Our calculations in § || show that accretion luminos- 
ity will emerge after recombination depletes the initial internal 
energy of the envelope, hence after a few tens of days. The 
corresponding luminosity at emergence will be ~ 10 37 ergs s , 
which with the present capabilities of HST would allow for de- 
tection at emergence perhaps out to 20-25 Mpc. The luminos- 
ity at emergence is higher than for SN1997D so that the galactic 
background would be less important, and ground based instru- 
ments might be relevant as well. 

How often can we expect an opportunity to observe a black 
hole emerging in a supernova? An upper limit is, of course, the 
rate at which black hole forming supernova e occur. Adopting 
the recent results of Fryer & Kalogera (1999) that about 10% of 
core collapse supernovae make black holes, we expect a rate of 
about one event per thousand years per galaxy. There are some- 
what less than one thousand galaxies at a distance <J 20 Mpc 
from the Milky way, and hence there may be approximately 
one observable event (equal or better than SN1997D) per sev- 
eral years. This estimate is an upper limit, since a significant 
fraction of high mass progenitors will end their lives as Type 
Ib/Ic supernovae after losing their hydrogen envelope, and such 
supernovae are unfavorable for a practical observation at emer- 
gence. Furthermore, even "perfect" candidates may not allow 
for an actual observation, as seems to have been the case in 
SN1994W. Following its recombination peak, the tail of the 
SN1994W light curve failed to settle on an exponential decay 
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that would arise due to 56 Co decay ( Sollerman et al., 1998 ). In- 
terestingly, the last data point led to a limit of < 2 x 10~ J on the 
mass of 56 Co in the ejected envelope. However, the light curve 
showed a very steep power-law in time until it dropped below 
the detection threshold some 220 days after the explosion. The 
luminosity during this phase was much too high to be du e to ac- 



cretion (10 40 - 10 42 ergs s" 1 ) and |Sollerman et al. (1998|) found 



that it was most likely caused by interaction with circumstellar 
material, with dust formation also playing an important role. 
We can only speculate that in the absence of circumstellar inter- 
action and dust extinction, the light curve would have declined 
exponentially after the recombination peak and the accretion 
luminosity may have emerged while still marginally detectable. 
The case of SN1994W implies that the rate of potential observ- 
able black hole emergence may be significantly smaller than 
that of black hole forming Type II supernovae. 

Nonetheless, we urge that any Type II supernovae which 



shows a diminished 56 Co abundance should be monitored, and 
its light curve be tracked for several months after the explosion. 
The prospects of observing black hole emergence may improve 
significantly with potential future instruments, such as NGST 
or a dedicated faint-supernovae project. 
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APPENDIX 

THE RESCALING SCHEME: RADIATION-HYDRODYNAMIC RENORMALIZATION 

The inclusion of the accreting region imposes an inherent obstacle in any numerical simulation of a supernova light curve. In the 
innermost accreting region the material is free-falling at a velocity which is a significant fraction of the speed of light, so that the 
dynamical time scale is less than 10~ 3 sec. On the other hand, the evolutionary time for the light curve is of the order of tens of days 
to clear the recombination phase, and years to reach black hole emergence if there is a non-negligible amount of radioactive elements 
in the envelope. Thi s extreme diversity of ti me scales prevents any single calculation which tracks the entire envelope. 



As described by Zampieri et al. (1998a), a considerable increase in the shortest physical time scale (and, hence, the Courant 
integration time step) can be obtained moving the inner boundary, /?,„, outward (at hundreds to thousands of Schwarzschild radii, but 
always well within the accreting region). This turns out to be possible because the gas in the inner accreting region is in LTE and in 
free-fall, thus allowing a reliable setting of the boundary conditions even if Rj„ is much larger than the black hole horizon. When a 
Lagrangian zone passes through the inner boundary, it is removed from the calculation and the remainder of the envelope quantities 
are rezoned to maintain a satisfactory resolution. 

Even with such a relocation of the inner boundary, the feasibility of a numerical computation requires some additional measures 
to accelerate the computation. To this end, the code incorporates the Multi-Time-Step Procedure, where the integration domain is 
separated into several subgrids, and each is evolved at its own local time scale; most of the CPU time is thus devoted to evolving 
the regions with smaller time steps. We note that this procedure requires a communication scheme between the sub-grids based on 
extrapolation, and can be safely used when there are no sharp features propagating through the envelope. Thus, we use it only after 
that the recombination front has completed its sweep of the envelope material. The speed-up allowed by the MTP is typically a factor 
of 5-8. 

Unfortunately, moving the inner boundary outward and using the MTP is not sufficient when the target physical time is of the order 
of years^J In the case of SN1997D, tracking the evolution of the envelope for ~ 1000 days would require roughly 1000 CPU hours 
on the SGI 500MHz EV6 Compaq XP1000 processor which we have used for the simulations. We were hence compelled to develop 
an additional acceleration algorithm, which is a rescaling scheme that allows us to integrate to the black hole emergence stage within 
a reasonable computational time. 

Principles of Rescaling 

The main objective in devising a rescaled model of a specific supernova evolution is to maintain all the quantitative features of the 
realistic model, while reducing the required computational time. The code uses an explicit Lagrangian scheme. Since the photon 
transport is treated by actually solving the transport equation, the reference velocity is the speed of light, and the Courant time step 
is determined by the width of the thinnest zone. This is always the innermost zone, AR,„(t ), which at any time t is a given fraction of 
the position of the inner boundary, Ri„(t). The efficiency of the code is essentially determined by comparing the (small) ratio of the 
time step to a characteristic evolutionary time - for example, the expansion time scale, to - 

AR in (t)/c R in {t)/c 

oc . (Al) 

to to 

The essence of a rescaling scheme would then be devising a rescaled model where 

Rl(t S )/c „ Rin(t)/c 

T s = a " » ( A2 ) 

9 We note that using a fully implicit code would not allow a significant gain in computational time because limitations on the time-step imposed by accuracy are 
comparable to limitations imposed by the Courant condition throughout most of the evolution. 
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where a subscript s denotes the rescaled quantities. The rescaled model is then more efficient by a factor of a > 1. 

Successful rescaling must (a) provide a simple relation between its own physical quantities and that of the realistic model and (b) 
maintain the hierarchy (inequalities) of time and energy scales that exists in the realistic model, so that the relative importance of the 
different physical process is carried on into the rescaled model. The rescaling scheme we have chosen to employ in our simulations 
is detailed below. 



The Rescaling Relations 

A natural rescaling scheme would be one where all the physical quantities of a mass element (r,v,T,p...) of the realistic model 
are mapped with simple linear relations to the rescaled model, i.e., r 1 = a B r, V s = ayv, etc. In our rescaling scheme we distinguish 
between the local thermodynamic and kinematic quantities of the mass elements in the envelope. Since the photon opacities are 
sensitive and complex functions of the thermodynamic quantities , i.e., density and temperature, we find it preferable to leave these 
quantities unchanged in the course of rescaling T s = T, p s = p, which leads to c s s = c s , k s = k. Note that such a choice also implies that 
the inner boundary of the realistic model should approximately map onto the inner boundary in the rescaled model, since the position 
of the inner boundary is determined by the requirement of LTE, which arises from the local values of temperature and density. 

Once this choice has been made, the constraints on the rescaling scheme are quite straight forward. If we require that the ratio of 
diffusion and recombination times scales to the expansion time scale remain unchanged, i.e, 

t to t Q tQ 

R s V s = RoVo. (A4) 



the rescaled model must satisfy (see eqs. iMpp): 



The condition ( A4) is sufficient to determine the rescaling strategy when consistently mapping all radii between the realistic and 
the rescaled models. If all radii are decreased by a factor of a, all velocities must be increased by the same factor; the efficiency of 

: a. 



the simulation (eq. [ A.1 )]) is then increased by a factor of a 1 /a 2 ■ 



If we further impose invariant luminosity ratios according to 

L diff,o _ L difffi L s rec 



L"acc,0 L acc fi L s accQ L aC c,0 

we arrive (using eqs. [EIR and Il9||) at the condition 



(A5) 



M s BH Rl=M BH Ro. (A6) 



The principle relations of our rescaling scheme are thus 



*oV'_ (m s bh 



RoJ \VoJ KM, 



(A7) 



This choice of rescaling allows us to consistently map the marginally bound shell of the envelope, since the marginally bound radius, 
R mb also rescales as R s mbfi = (2M BH (t^ l ) 2 ) 1 ^ = a _1 ^ mi ,o. 
Further relations that immediately follow are 



-,)-" •\jZ)=° '\j£)- a •\sr a ■ <A8) 

The efficiency of rescaling is explicit in the rescaling of time - in the rescaled model all physical times are reduced by a factor of 
a 2 , while the numerical time step is only reduced by a factor of R s in /Rj„ = a. 

Limitations of Rescaling 

The rescaling scheme presented above is not without limitations; most notably, the acceleration factor, a is limited to a value 
of a few, as discussed below. Nonetheless, even a value of a = 5, which was used in the case of SN1997D reduces the required 
computational time from weeks to days, hence allowing for a practical numerical investigation. 



The Acceleration Factor 

Since the rescaling scheme calls for an increase of a factor of a of the envelope velocities, clearly a cannot be arbitrarily large, 
since we must maintain a sub-relativistic envelope. Furthermore, since the black hole mass is increased, its Schwarzschild radius 
increases as well, and naturally it should not be allowed to approach the integration domain throughout the calculation. Both these 
constraints result in the practical limit on the value of a to no more than a few. 
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Radioactive Elements 



The scaling laws of luminosity and mass in equation ( A8) create an undesired inconsistency regarding the luminosity from ra- 
dioactive heating, since for this source LocM. In order to recover the correct rescaled power output from radioactive decays, the 
energy generation rate from radioactive decays per unit mass in the envelope must therefore be increase by a factor of a 2 . This 
could potentially impose a significant deviation in the evolution of the rescaled model from that of the original one, if it artificially 
enhances the importance of energy deposition from radioactive decays when compared to the other energy scales. However, since 
we limit ourselves to envelopes where the abundance of radioactive elements is very small, the artificial amplification of the energy 
production per unit mass of radioactive decays has a negligible effect. 

It is further noteworthy that the effective transparency of the photons to the 7-ray photons emitted in the decays must be adjusted 
in order to recover the appropriate scaling of the 7-transparency time. The simple form of 7-ray opacity (eq. [lVA) allows to achieve 
consistency with the change = q;k 7 , so that tf rans /?q = t tmns ^/tQ. A similar transformation for the optical photons would be more 
complex, but can be neglected since the onset of thermal-photon transparency is imposed by recombination (rather than expansion), 
and is chiefly determined by the thermodynamic variables p and T, which are kept invariant. 

The Accretion Time Scale and Radius 

A second inconsistency that arises in this particular rescaling scheme concerns the initial accretion time and radius, t acc o, R a c C ,o- As 
is evident from equations ( |l3| , [T^ ), when the thermodynamic properties of the helium layer are not changed, so that the sound speed 
is conserved, t acc ,R acc oc Mbh- Our rescaling scheme actually increases these quantities by a factor of a, while other characteristic 
times and radii are decreased by factors of a 2 and a, respectively. This is an inevitable consequence of our choice of rescaling. 

This inc onsistency is significant, however, only while the accretion flow is Bondi ( 1952)-like, when t acc < tp and R acc < R,„b (C olpi 
et al., 1996). As the envelope continues to expand, the hierarchy of time scales and radii is eventually reversed and the accretion 
eventually transforms into a dust-like flow. Since the accretion time does not enter explicitly in any of the estimates for evolution 
and luminosity (see § |J), it is only this hierarchy which is of real consequence, and hence rescaling can be applied safely to models 
which are characterized by t acc ^> to, R aC c ^> Rmb- 

Kinetic vs. Thermal Energies 

A similar analysis concerns the ratio of kinetic and thermal energies in the envelope. It is clear that the kinetic energy will scale 
as Ekin ~MV 2 — + E-l jn /Eki„ = a~ x , while the internal energy (neglecting radioactive heating) will scale as E tn ^R 3 — > EfjE t h = a~ 3 ', 
since the temperatures are unchanged. Rescaling thus amplifies the relative importance of kinetic energy over thermal energy in the 
expanding envelope. 

However, as is the case for the accretion and expansion time scales, the real issue is the hierarchy of energy scales. In fact, 
the assumption of homologous expansion is equivalent to assumption that the conversion of thermal energy to kinetic energy is 
negligible, i.e., that Ekin ^ E tn , which is satisfied in the context of our rescaled models. 

A Test of Rescaling 

We demonstrate here the reliability of our rescaling scheme by studying a specific example. We assume a supernova envelope 
somewhat similar to the SN1997D case examined in § |], with sev eral simp lifications that allow for a more rapid computation. The 



basic features of the model are presented in table A2 and in Figs. A.4 - A.4. We note that the initial velocity profile throughout the 
entire envelope was set so that v(r) = r/to. The initial accretion time-scale is determined by the properties of the region of bound 
material. A single fiducial radioactive isotope was included, with a mass fraction which is 0.01 of that of oxygen. This isotope was 
assumed to emit all the decay energy in 7— rays. 
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FIG. All. — Initial profile for the toy model: density, tem- 
perature and velocity. 



We have carried out two simulations of this model, the first in its original form and once rescaled by a factor a = 5 (sin ce fp < C t acc fi 
and Rmbfl < Raccfl to begin with, the model can be rescaled at t = 0). The light curves of both models are compared in Fig. A13| , where 
luminosity and time in the rescaled model were rescaled inversely in order to compare them to the results of the original model, i.e., 
L s (t s ) — > aL s (a 2 t s ), . We note that the rescaled model required only 60 CPU hours, compared to about 300 CPU hours for the original 
model to reach an equivalent evolutionary stage. 



Table A2 

Parameters for the toy model used to test rescaling * 



model 


Mbh 


M M 


U B ,' 




Rom 




V au ,' 






r ra h 




to 


'accfi 




(Me) 


CM©) 


(Af ) 


(Mo) 


(cm) 


(cm) 


(cms 4 ) 


(Mq) 


(ergs gm _1 s~') 


(days) 


(cm 2 gm _1 ) 


(hrs) 


(hrs) 


Original 


3 


18.75 


3.75 


15 


2 x 10 13 


4.5 x 1() 12 


6.67 x 10 s 


0.036 


9 X 10 s 


50 


0.01 


25/3 


50 


Rescaled 


15 


0.15 


0.03 


0.12 


4 x 10 12 


9 x 10" 


3.35 x 10 9 


2.88 x 10~ 4 


2.25 x 10 10 


2 


0.05 


1/3 


250 



*Rescaled by a factor a = 5 (see text) 
a mass of helium-rich layer 
b mass of hydrogen rich layer 
c outer radius 

d outer radius of helium-rich layer 

e initial velocity at outer radius 

f total mass of radioactive isotope in the envelope 

g energy generation rate per unit mass of radioactive isotope 

h life-time of radioactive isotope 

'7-ray opacity per electron radioactive isotope 
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The reliability of rescaling is evident from the very good agreement of the light curve and accretion rate history. In particular, the 
accretion luminosities at emergence agree to better than 10% (and a similar agreement is found in the accretion rates). We note that 
there is an obvious deviation in the light curve during the recombination phase, which is not unexpected since recombination is a 
non-linear process ( \rnett, 1996): the rescaled model has smaller physical size, and the recombination front crosses it more quickly. 
This discrepancy is not important for determining the late-time accretion rate and luminosity, since the accreting region is unaffected 
by the recombination occurring much further out in the expanding envelope. 
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fcfj 
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original model 
rescaled model 




log 10 t (days) 

FIG . A 1 3 . — Calculated bolometric light curve of the toy model shown in Fig. [a4| (solid lines), compared with the calculated light curve of model rescaled by a 
factor of 5 (dashed lines). The inset shows the light curves at the time of emergence of the black hole 
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